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LONG TIME DYNAMICS AND DISORDER-INDUCED TRAVELING 
WAVES IN THE STOCHASTIC KURAMOTO MODEL 

ERIC LUgON AND CHRISTOPHE POQUET 


Abstract. The aim of the paper is to address the long time behavior of the Kuramoto 
model of mean-field coupled phase rotators, subject to white noise and quenched frequen¬ 
cies. We analyse the influence of the fluctuations of both thermal noise and frequencies 
(seen as a disorder) on a large but finite population of N rotators, in the case where the 
law of the disorder is symmetric. On a finite time scale [0, T], the system is known to be 
self-averaging: the empirical measure of the system converges as N —► oo to the deter¬ 
ministic solution of a nonlinear Fokker-Planck equation which exhibits a stable manifold 
of synchronized stationary profiles for large interaction. On longer time scales, compe¬ 
tition between the finite-size effects of the noise and disorder makes the system deviate 
from this mean-field behavior. In the main result of the paper we show that on a time 
scale of order V7V the fluctuations of the disorder prevail over the fluctuations of the 
noise: we establish the existence of disorder- induced traveling waves for the empirical 
measure along the stationary manifold. This result is proved for fixed realizations of 
the disorder and emphasis is put on the influence of the asymmetry of these quenched 
frequencies on the direction and speed of rotation of the system. Asymptotics on the 
drift are provided in the limit of small disorder. 


1. Introduction 

1.1. Long time dynamics of mean-field interacting particle systems. The macro¬ 
scopic behavior of numerous stochastic interacting particle systems appearing in physics 
or biology is usually described by nonlinear partial differential equations. In this context, 
systems of diffusions in all-to-all interactions, that is mean-field particle systems [32} [33], 
have attracted much attention in the past years, since they are relevant in many situations 
from statistical physics (synchronization of oscillators [1} [271 , 4Tj) to biology (emergence 
of synchrony in neural networks M) and have provided particle approximations for var¬ 
ious PDEs (see nn no] and references therein). From a statistical physics point of view, 
a natural extension of these models concerns similar particle systems in a random envi¬ 
ronment, that is when the particles obey to the influence of an additional randomness, 
or disorder, representing inhomogeneous behaviors between particles. Such a modeling 
is particularly relevant in a biological context, where each particle/diffusion captures the 
state of one single individual (activity of a neuron, phase in a circadian rythm) and the dis¬ 
order models intrinsic dynamical behavior for each individual (e.g. inhibition or excitation 
in populations of heterogeneous neurons M)- 

The aim of the paper is to address the influence of the disorder on the long time dynamics 
of a large but finite population of mean-field interacting diffusions with noise. A crucial 
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aspect in this perspective is the notion of self-averaging: in the limit of a large number of 
individuals and/or on a long time scale (in a way that needs to be made precise), is the 
macroscopic behavior of the system the same for every typical realization of the disorder? 
If not, is it possible to quantify the influence of the fluctuations of the random environment 
on the behavior of the system? 

It appears that the analysis of such mean-held systems differs significantly depending on 
the time scale one considers. On a time scale of order 1 (w.r.t. the size of the population), 
it is now well-known that the macroscopic behavior of mean-held particle systems are well 
described by nonlinear PDEs of McKean-Vlasov type [201 M- A vast literature exists 
on the links between the microscopic system and its mean-held limit (fluctuations, large 
deviations and finite time dynamics) mostly in the non-disordered case (see e.g. HhI-TT . 13 
and references therein) but also for disordered systems [161 [28] . 

When one considers longer time scales (w.r.t. the size of the population) and for a large 
but finite number of particles, some randomness remains in the system so that Brownian 
fluctuations generally induce microscopic dynamics that may differ signihcantly from the 
dynamics of the mean-held equation. For mean-held systems without disorder, a vast 
literature exists concerning huctuations induced by thermal noise. In this respect, the 
notion of uniform propagation of chaos has been addressed for several mean-held models 
by many authors (see e.g. EHi for the granular media equation or [25 1 [39] for ranked- 
based models). In case the mean-held PDE admits an isolated stable hxed point, due to 
large deviation phenomena, the finite-size system exits from any neighborhood of the hxed 
point at exponential times in N ( N being the size of the population) [171 35] . whereas 
in case of an unstable hxed point, the system escapes at a time scale of order log N [38] . 
Fewer results exist in the case where the mean-held PDE admits a whole stable curve of 
stationary solutions. In [?, T5j, the effect of thermal noise is considered for the mean-held 
plane rotators model [6j which is known to admit in the limit asJV->c»a stable circle of 
stationary solutions. In this case, the hnite size particle system has Brownian huctuations 
on time scales of order N. 

In the case of disordered systems, we are not aware of any similar analysis on long 
time dynamics of mean-held interacting particles. The present work could be seen as a 
first result in this direction. In particular, we provide in Theorem 12.31 a rigorous and 
quantitative justihcation to a phenomenon already observed by Balmforth and Sassi [2] 
on the basis of numerical simulations. 

1.2. The stochastic Kuramoto model with disorder. We address in this paper the 
long time behavior of the Kuramoto model with noise and disorder, which describes the 
evolution of a population of rotators (the j th rotator being dehned by its phase gV (t) € 
T := M/27tZ), given by the system of N ^ 1 stochastic differential equations of mean-held 
type 

I\ 

d <Pj(t) = Sujj dt - (*) “ (^f(t))dt + adBj{t), j = 1,...,N, t > 0, (1.1) 

v i=i 

where n is a family of standard independent Brownian motions, K , a and 5 are 

positive parameters. In particular, <5 > 0 is a scaling parameter. The main result will be 
stated for small 5 > 0, as it relies on perturbation results of the case where <5 = 0. 

The Kuramoto model niETiiii] is the main prototype for synchronization phenomena 
and, due to its mathematical tractability, has been studied in details in the past years 

mmmm- 
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Remark 1.1. Note that (EH) is invariant by rotation: if (i))j=i,...,Ar solves ED. then 
so does {gj^if) + a)j=i,...,N for all a € M. Moreover, by the change of variables t — > t/a 2 , 
one can get rid of the coefficient a in front of the Brownian motions (up to the obvious 
modifications 5 —>• 5/a 2 and K —>• K/a 2 ). Hence, with no loss of generality, we suppose 
a = 1 in the following. 

Following the point of view adopted at the beginning of this introduction, the system 
ED presents two types of noise: in addition to the thermal noise (Bj), the disorder in 
ED is given by a sequence (uj)j=i,...,N of i.i.d random variables with distribution A, 
independent from the Brownian motions. Each ojj represents an intrinsic inhomogeneous 
frequency for the rotator gfij. The index uj in the notation gj^ is used to emphasize the 
dependency of the system in the disorder. 

A crucial aspect in the understanding of the dynamics of ED concerns the (possible 
lack of) symmetry of the sequence (ujf)j ^ \. First note that, by the obvious change 
of variables gfij(t) i-a <g>f (t) — E(w)f in (II.ID . it is always possible to assume that the 
expectation of the disorder E(w) = J R o;A(da;) is zero (otherwise, we observe macroscopic 
traveling waves with speed E(cn)). The asymmetry of the disorder can be given at different 
scales. The most simple situation corresponds to a macroscopic asymmetry, that is when 
the law A itself is asymmetric. With no loss of generality, we can for example assume that, 
on a macroscopic level, a majority of rotators will be associated to a positive frequency 
whereas a minority will have negative frequencies. In the limit of an infinite population, 
this asymmetry makes the whole system rotate at a constant speed that depends only 
on the law A and this rotation is noticeable at the scale of the nonlinear Fokker-Planck 
equation (I1.3P associated to (11.11) . This case has been the object of a previous paper (see 
m , Theorem 2.2 and Section [2.21 below). 

The present paper is concerned with the situation where the law of the disorder is 
symmetric. Here, the previous argument cannot be applied since in the limit as N —>• 
oo, the population is equally balanced between positive and negative frequencies: the 
macroscopic speed of rotation found in [21], Theorem 2.2 vanishes. Hence, the analysis of 
long time dynamics of ED requires a deeper understanding of the microscopic asymmetry 
of the disorder, that is the finite-size fluctuations of the disorder w.r.t. the thermal noise. 
An informal description of the dynamics of ED is the following (see Figure [2] below): if the 
constant K is sufficiently large, the mean-held coupling term leads to synchronization of 
the whole system along a nontrivial density. Even if A is symmetric, finite-size fluctuations 
of the sample (wj)j=i,...,jv make it not symmetric so that the fluctuations of the disorder 
compete with the fluctuations of the Brownian motions (Bj)j = i n and make the whole 
system rotate with speed and direction depending on the fixed realization of the disorder 
(ojj) (and not only on the law A itself). The main point of the paper is to give a rigorous 
meaning to this phenomenon, noticed numerically in [3]: we will show that at times of 
order VlV, the dynamics of (11.11) deviates from its mean-held limit, with the apparition 
of synchronized traveling waves induced by the finite-size huctuations of the disorder. We 
refer to Paragraph 11.61 below for a precise description of this phenomenon. 

We present in the following subsections some well-known properties of (11.11) which are 
needed to state our result. We describe in particular its infinite population limit on 
bounded time intervals and the existence of stationary measures for the limit system in 
case of symmetric disorder. 

1.3. Mean-field limit on bounded time intervals. All the statistical information of 
ED is contained in the empirical measure (n% t )t > o € C([0,oo), Aii(T x M)) (M.\ being 
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the set of probability measures endowed with its weak topology) defined as 

1 N 

VN,t Jt t^O. ( 1 . 2 ) 

3 = 1 

When the distribution A of the disorder satisfies /|u;|A(du;) < oo and the initial con¬ 
dition /i^ 0 converges weakly to some po when N —» oo, it is easy to see ([16, [28] ) 
that the empirical measure (11 .2 jl converges weakly on bounded time intervals (that is 
in C'([0, T], A4i(T x M)) for all T ^ 0) to a deterministic limit measure whose density pt 
with respect t <8> A (where l denotes the Lebesgue measure on T) satisfies the following 
system of nonlinear Fokker-Planck PDEs: 

d t p t (9,uj ) = ^d%p t (9,u) - d e (p t (9,uj)((J *p t ) x (9) + Su)^, u € Supp(A), 6 € T, t>0, 

(1.3) 

where 

J(0) : = - Ksm(9 ), (1.4) 

and (-)a represents the integration with respect to A: (J * u)\(6) = f R f T J(^)u(9 — 
dipX(du). We insist on the fact that in (11.31) . a; is a real number in the support of 
A, while in m and C2D, it is an index emphasizing the dependency in the disorder of 
the system. 

Some properties of system (|1.3I) are detailed in [21]. In particular, if A-almost surely, 
Po(-,uj) is a probability measure then (|1.3I) admits a unique solution p t for all t > 0 such 
that A-almost surely, pt(-,uj) is also a probability measure, with positive density with 
respect to the Lebesgue measure and is an element of C'°°((0,oo) x T, M). 

1.4. Symmetric disorder. As already mentioned, we consider the case where the law 
A of the disorder is symmetric. We restrict our analysis to finite disorder: fix d ^ 1 and 
suppose that the frequencies (ojj)j ^ i take their values in {u~ d , ..., uj d }, 

where a;* = — w - * for all i = 0, We denote as (A* G [0,1],* = —d,...,d) the 

probability of drawing each c o l and assume that A* = A - * for all * = 1,..., d. From now 
on, the law of the disorder A is identified with (A~^,..., A 01 ). Note that we may suppose 
in the following that ujq = 0 ^ Supp(A). The result still holds with obvious changes in 
notations. 

Under this hypothesis, almost surely, for sufficiently large N, each possible value w* of 
the disorder appears at least once and we can rewrite m by regrouping the rotators 
into (2d + 1) sub-populations: for all i = —d ,..., d, denote as N z the number of rota¬ 
tors (ipj with frequency uj 1 . Obviously, N = Yli=-d,N l and the system (11.11) 

becomes 

d N k 

d <p)(t) = 5u l dt-^ ^sin(^(t)-^f(t))dt+dS](t), j = l,...,iV\ i = -d,...,d. 

k=—d 1=1 

(1.5) 

In this framework, the empirical measure t in (11.21) can be identified with (p N d t , • • •, t ), 
where p l N is the empirical measure of the rotators with frequency oj 1 : 
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and its mean-field limit (11.31) can be identified with pt = (p t d ,... ,pf ), solution to 

dtpW) = ld 2 0 pi(e)-de^pi(9)^ X k J*p^e) + Su J ^, o, i = -d,...,d. (1.7) 

1.5. Stationary solutions and phase transition. A remarkable aspect of the Ku- 

ramoto model is that one can compute semi-explicitly the stationary solutions of (II.7j) . 
when A is symmetric (see e.g. my each stationary solution to ([1.711 is the rotation of a 
profile q = (q ~ d ,..., q d ) (i.e. given by q(- + a) for some a € T)) of the form 

S\{e,2Kr) 


q\6) = 


where for each i = —d,. 


Z' s (2Kr) 

,d, q*(-) is a probability density on T. S$(6,2Kr) is given by 


( 1 . 8 ) 


SK6,x) = 


_ xcos9+25u) l 6 


(1-e 


r6 

Jo 


—x cos u—25u)’ L z 


d u 


+ e 


AnSuj 1 


r 2?r 


—x cos u—2buS l u 


du 


(1.9) 


?!(2 Kr) is a normalization constant and r is a solution of the hxed-point problem 


r = d' < s(2 Kr) 


with 


*«(*) = ^ A* 


/o cos(0)Sg(0,x) d6 


Z\ 


(*) 


( 1 . 10 ) 


( 1 . 11 ) 


k=—d 

We refer to m or [25], p. 75 for more details on this calculation. Computing the so¬ 
lution to the fixed-point relation (11.101) enables to exhibit a phase transition for (11.71) : 
the value r = 0 always solves (11.101) and corresponds to the uniform stationary solution 
q = (1/27r,... , l/27r). It is the only stationary solution to (11.71) as long as K ^ K c , for a 
certain critical parameter K c = K c (5, (w*)*, (A*),) > 1. This characterizes the absence of 
synchrony in case of small interaction. When K > K c , this flat profile coexists with circles 
of synchronized solutions corresponding to positive fixed-points in (jl.lOD : each solution 
r > 0 to (jl.101) gives rise to a nontrivial stationary profile q given by (|1.8|) and to the circle 
of all its translation q(- + a), by invariance by rotation of the system (see Figured]). 

However, several circles may coexist when K > K c and these circles may not be locally 
stable (even the characterization of these circles in full generality is unclear, see e.g. [29], 
§ 2.2.2). To ensure uniqueness and stability of a circle of non-trivial profiles, fix K > 1 
and restrict to small values of 5: it is indeed proved in m , Lemma 2.3 that there exists 
<51 = 5\{K) > 0 such that for all 5 ^ <5i, the fixed-point problem (11.101) admits a unique 
positive solution r$. We denote by the corresponding profile given by (11.81) with r = rg, 
by q^p t s its rotation of angle ij> £ T (i.e. 5^5(•) := qo t s(- — i/O) and by M the corresponding 
circle of stationary profiles (see Figured]): 

M := {q^s : iff € T} . (1.12) 

It is proved in |2T|, Theorems 2.2 and 2.5 that the circle M is stable under the evolu¬ 
tion (ED: the solution of O) starting from an initial condition sufficiently close to M 
converges to a element of M as t —>• 00. More details about this stability are given 
in Section 12.31 Whenever it is clear from the context, we will use the notations qs or qy, 
instead of q^ t s, depending on the parameter we want to emphasize. 
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(a) Correspondance between 
fixed-points of $a(-) and sta¬ 
tionary solutions to GZ3- 


(b) A synchronized profile with 
d=2,q = (q~ 2 , 9 _1 , q 1 , q 2 )- 


Figure 1. Fixed-point function and stationary profiles when K = 5, d = 2, coi = 1, 
102 = 10 and 5 = 0.1. 
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Figure 2. Evolution of the marginal of the empirical measure (11.211 on T for a fixed 
choice of the disorder (N = 600, A = l(<5_i + <5i), I\ = 6). Starting from uniformly 
distributed rotators on T (t = 0), the empirical measure converges to a synchronized 
profile on the manifold M (t = 6) and then moves (here to the right) at a constant 
speed, on a time scale compatible with N 1 ^ 2 . 

1.6. Long time behavior. Simulations of (11.51) (Figure [2]) suggest an initial transition 
of the system from an incoherent state to a synchronized one, during which the empiri¬ 
cal measures of the rotators approaches the circle M of synchronized stationary profiles. 
Secondly, the empirical measure remains close to M and travels at first order at constant 
speed (which is random, depending on the realization of the disorder, see Figure [3]) along 
M on the time scale lV 1//2 t. Let us give some intuition of this phenomenon: to fix ideas, 
consider the case where d = 1, uj 1 = — = 1 and A -1 = A 1 = ^. This corresponds to 
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Figure 3. Trajectories of the center of synchronization for different realizations of the 
disorder (A = i(5_o.s + <5o.s), K = 4, N = 400). 


the simplest decomposition in (11.511 between two subpopulations, one naturally rotating 
clockwise (c Oi = +1) and the second rotating anti-clockwise (c Oi = — 1). One can imagine 
that fluctuations in the finite sample (uq,..., uqv) € {±1}^ may lead, for example, to a 
majority of +1 with respect to —1, so that the rotators with positive frequency induce a 
global rotation of the whole system in the direction of the majority. When N is large, this 
asymmetry is small, typically of order 1W 1 / 2 and is not sufficient to make the empirical 
measure drift away from the attracting manifold M, but induces a small drift that becomes 
macroscopic at times of order IV 1 / 2 . 

The purpose of the paper is precisely to prove the existence of this random traveling 
wave and show that it is indeed an effect of the fluctuations of the disorder. Our approach 
consists in a precise analysis of the dynamics of the empirical measure (11.61) , which involves 
both disorder and thermal noise. One of the main difficulties is to control the thermal 
noise term and prove that it does not play any role at first order on the A rl / 2 -time scale. 

2. Main results and strategy of proof 


2.1. The result. 


Admissible sequence of disorder. We stress the fact that the random traveling waves de¬ 
scribed above is essentially a quenched phenomenon, that is, true for a fixed realization of 
the disorder (uii)i ^ i. In particular, the result does not really depend on the underlying 
mechanism that produced the sequence {uii)i ^ i, it only depends on the asymmetry of this 
sequence. We prove our result for any admissible sequence of disorder {uf)i ^ i, defined as 
follows. 


Definition 2.1. Fix a sequence (aq)j ^ i taking values in {uj d ,ui *- d l \ ... ,uj d 1 ,co d } and 
for all N ^ 1, define the empirical proportions of frequencies in the N-sample (uq,... ,un) 






( 2 . 1 ) 
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where N k is the number of rotators with frequencies equal to uj k (recall Section \l.f\ ). Define 
also the fluctuation process associated to (ad ^ i by := ■ ■ ■, where 

d := JV 1/2 (A& - A fc ), k = -d,...,d, Nfil, (2.2) 

where ( X ~ d ,..., X d ) is given in Section [TfJ j Note that Ylk=-d^N = 0 f or all N fi 1. We 
say that the sequence (a;*),; ^ i is admissible if the following holds 

(1) Law of large numbers: for all k = —d ,..., d, X^ converges to X k , as N —» oo. 

( 2 ) Central limit behavior: for all f > 0, there exists Nq (possibly depending on the 
sequence (> i) such that for all N fi Nq, 


max 

k=—d,...,d 



^ 7V C . 


Remark 2.2 (Admissibility for i.i.d. variables). An easy application of the Borel-Cantelli 
Lemma shows that any independent and identically distributed sequence of disorder (u ;$)$ j> i 
with law X is almost surely admissible, in the sense of Definition \2. A 


Main result. From now on, we fix once and for all an admissible sequence (wj)j ^ i in the 
sense of Definition 12.11 A convenient framework for the analysis of the dynamics of (|1.6I) 
and (HAD corresponds to the space H d 1 , dual of the space H d , which is the closure of the 
set of vectors ( u ~ d ,..., u d ) of regular functions u k with zero mean value on T under the 
norm 


1,<2 : — 



{d e u\e)) 2 do 


1/2 


(2.3) 


Remark that if u is a vector of probability measures on T, then u naturally belongs 
to Hf l , since the family of vectors given by a n> fc( 0 ) = ( 0 ,..., 0 , , 0 ,..., 0 ) and 

b n ,k(8) = (0,... , 0, , 0,... , 0) form an orthonormal basis of H d and for each such 

vector u 


N -i,d = 


\ 


y y ({u, a n ,k ) 2 + {u, dfc) 2 ) ^ 7T 


k=—dn= 1 


\ 


* E a k r' 


(2.4) 


k=—d 


More details on the construction of H d 1 are given in Appendix [A] The main result of the 
paper is the following. 


Theorem 2.3. For all K > 1, there exists 5(K) such that, for all 5 ^ S(K), there exists 
a linear form b : R M+1 —>• R (depending in K, 5, the probability distribution X and the 
possible values of the disorder u i l ) and a real number £q > 0 such that, for any admissible 
sequence (cjj)i y i, any vector of probability measures po satisfying dist H -i(po, M) ^ £o 
such that for all e > 0 , 

p (\\hN ,0 Po||_! ^ > e) —t 0, as N -A oo , (2.5) 

then, there exists 9 q £ T (depending on po) and a constant c such that for each finite time 
tf > 0 and all e > 0, denoting t{f = cN logiV, we have 


sup 

■o ■*/) 


K te[tfi 


~ fte 0 +b(t N )t 


-1 ,d 


fie) —> 0, as N —>• oo . 


( 2 . 6 ) 
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Moreover, £ i-a &(£) has the following expansion in 6: for all £ such that Ylt=-d^ k = 
we have 

d 

6(0 = -5 ]T ^ k + 0(5 2 ). (2.7) 

k=—d 

Theorem 12.31 is simply saying that, on a time scale of order N 1//2 , the empirical measure 
m is asymptotically close to a synchronized profile q € M, traveling at speed &(£jvO 
along M. This drift depends on the asymmetry of the quenched disorder (wj)i^i. 
In (12.61) . t(f represents the time necessary for the system to get sufficiently close to the 
manifold M. 

Some particular cases and extensions. First remark that the situation where the sample of 
the disorder ( 0 ^)*=!....^ is perfectly symmetric corresponds to for alH = 1,..., d. 

In this case, the drift in (12.61) vanishes: 

Proposition 2.4. If for all i = 1,..., d, , then b{ff) = 0. 

In particular, if one chooses the disorder in such a way that is always symmetric 

(e.g. choose an even number of particles N and define each Wj to be alternatively ±1), 
the drift is always zero. We believe in this case that one would need to look at larger 
time scales of order N to see the first order of the expansion of the empirical measure pjsr. 
Proof of Proposition 12.41 is given in Section 17.11 

In case the sequence (ojfli 1 is i.i.d. with law A, a standard Central Limit Theorem 
shows that the drift 6 (£at) converges in law to a Gaussian distribution Af(0,v 2 ), where v 2 
depends on K, 6, the probability distribution A and the possible values of the disorder uv l . 

Proposition 2.5. The following asymptotic of v 2 holds when 6 —>• 0: 

d 

v 2 = 5 A fc (c v k ) 2 + 0(d 2 ). (2.8) 

k=—d 

Proof of Proposition 12.51 is given in Section 17.21 

Remark 2.6. Without much modification in the proof, the result can be easily extended 
to sequences ( ujf)i ^ 1 with fluctuations of order different from V~N, that is when for some 
a € (0,1), 

—>• £ a , as N —>• 00 , (2.9) 

for some vector £ a where fffl := N a { Aat —A). In this case, the correct time renormalization 
is N a and we obtain a result of the type 

P sup \\pN,N°-t - g 0 o + 6 (f“)t||_i d > e \ -4-0, as IV 00 . (2.10) 

4/1 / 

Here, we only treat the case a = 1/2 for simplicity. For smaller fluctuations of size N~ a 
with a ^ 1, the time renormalization should be of order N. Since at this scale the effects 
of the thermal noise appear, the limit phase dynamics should be of diffusive type and a 
precise analysis of the different terms and symmetries that occur would be necessary to get 
the proper drift in this case. 


2.2. Links with existing models. 



10 


ERIC LUgON AND CHRISTOPHE POQUET 


Symmetric versus non-symmetric disorder. This work is the natural continuation of m, 
Theorems 2.2 and 2.5 in the case of a symmetric disorder. The purpose of [21] was to 
analyze the dynamics of the nonlinear Fokker-Planck equation (11,31) for both symmetric 
and asymmetric law of the disorder. The main point is that understanding (jl.3l) is not 
sufficient in itself for the analysis of the finite size system m in the symmetric case, 
since it does not account for the finite-size effects of the disorder that are crucial here. 

As already mentioned, in the case where A is asymmetric, one observes macroscopic 
travelling waves with deterministic drift at the scale of the nonlinear Fokker-Planck equa¬ 
tion (11.31) . It is reasonable to think that an analysis similar to what has been done in this 
paper would also show the existence of a finite order correction to this deterministic drift 
for a large but finite system with quenched disorder. 

Some previous results already suggested the possibility of these disorder-induced trav¬ 
eling waves in the Kuramoto model. Namely, the purpose of previous work [28] was to 
prove a quenched fluctuation result for the empirical measure (11.21) around its mean-field 
limit (| 1.3 D on a finite time horizon [0,T]. The main conclusion of [28] was that these 
fluctuations are disorder dependent and the long time analysis of the limiting fluctuations 
m suggested a non-self-averaging phenomenon for m similar to the one observed here. 

The case <5 = 0. This paper uses techniques previously developed in p 7 ] in the context of 
the stochastic Kuramoto model without disorder, that is when one takes 5 = 0 in (11.11) : 

K 

d <Pj(t) = -—^2sin(ip j (t)-ipi{t))dt+ dBj(t), j = l,...,N, (2.11) 

v i=i 

associated in the limit N — > oo to the mean-field PDE 

dtPt{0) = ^dgp t {6) - d e (pt(0)J *Pt(d)') ■ (2.12) 

Similarly to (11.71) in Section [L5l evolution ()2.121) generates a stable circle M 0 of stationary 
synchronized profiles when K > A' c (0) = 1 (see Section [BT] for further details). The model 
(12. 11 D - (12. 121) has been the subject of a series of recent papers [6, 0, .22, 23], addressing the 
linear and nonlinear stability of the circle of synchronized profiles Mq as well as the long 
time dynamics of the microscopic system (12.111) . The analysis of (12. 12|) strongly relies on 
the reversibility of (12.111) (with the existence of a proper Lyapunov functional, see [6] for 
more details), whereas reversibility is lost when 6 > 0. 

Concerning the long time behavior of (12.111) . it is shown in [TJ that under very general 
hypotheses on the initial condition, the empirical measure of (12.111) first approaches the 
circle Mq exponentially fast (that corresponds to the synchronization of the system (12.111) 
along a stationary profile solving (|2.121) 1 and then stays close to Mq for a long time with 
high probability, while the phase of its projection on Mq performs a Brownian motion as 
N —>• oo which corresponds to a macroscopic effect of the thermal noise. The persistence of 
proximity of the empirical measure to M for long times and the convergence of this phase 
to a Brownian motion were in fact already established in the unpublished PhD Thesis m 
the authors of (7J were not aware of, using in particular moderate deviations estimates of 
the mean field process. Note that the techniques of m do not apply here, since a similar 
analysis would involve moderate (or large) deviations in a quenched set-up, result that, to 
the best of our knowledge, has not been proven so far (for averaged large deviations, see 

A significant difference between [71 [15] and the present analysis is that the Brownian 
excursions in mm occur on a time scale of order N whereas it is sufficient to look at times 
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of order N 1//2 to see the traveling waves in the disordered case. This will entail significant 
simplifications in the analysis of (|1.5|) . since the detailed analysis on the thermal noise 
performed in mm will not be required here. 

Note also that, contrary to [7], we do not prove the first step of the phenomenon 
described in Figure El that is the initial approach of the system to a neighborhood of 
the manifold M in an exponentially short time, regardless of the initial condition. This 
result would require a global stability result for the system of PDEs (II.7[) which has not 
been proved for the moment, due to the absence of any Lyapunov functional for (11.71) 
when 5 > 0. We prove our result for initial conditions belonging to some macroscopic 
neighborhood of M (see Section [6] for more details). 

SPDE models with vanishing noise. This paper is related to previous works in the context 
of SPDE models for phase separation. In [12, [T9], the authors studied the Allen-Cahn 
model with symmetric bistable potential and vanishing noise. They showed that for an 
initial data close a profile connecting the two phase, the interface performs a Brownian 
motion. Some techniques initially introduced in these works, as the discretization of the 
dynamics in an iterative scheme, were developed in [7] in the context of the Kuramoto 
model without disorder (making use of Sobolev spaces with negative exponents to deal 
with empirical measures) and will have a central role in our analysis (see Section ESI). 
The results of [121 19] have been extended in [IT] by considering small asymmetries in the 
potential which induce a drift in the interface dynamics and by considering macroscopi- 
cally finite volumes [5j, with effect a repulsion at the boundary for the phase. Stochastic 
interface motions have also been recently studied in the context of the Cahn-Hilliard model 
with vanishing colored noise |2J. In this model, the limit behavior of the interface is given 
by a SDE (or system of SDE’s in the case of several interfaces) with drift and diffusion 
coefficients depending on coloration of the noise and on the length of the interface. 


2.3. Linear stability of stationary solutions. In the whole paper, we suppose that 
K > 1 and that <5 > 0 is smaller than some 5(K ) > 0. This critical value 5(K) is 
determined by S(K) = min(<5i(A'),^(it)), where <5i(/i) ensures the existence of a unique 
circle M of stationary solutions (recall Section II.5[) and where S- 2 (K) comes from the 
stability analysis of this circle (see Appendix [B] for more details). 

More precisely, our result relies deeply on the linear stability of the dynamical system 
induced by the limit system of PDEs (11.71) in the neighborhood of the circle of stationary 
profiles M. For ip € T, <5 > 0, consider the operator of the linearized evolution around 
qp,s € M given by 


= ^dgu 1 - 5uj l d e u l - d e 
for all i = — d ,..., d with domain 

ju = (u~ d ,...,u d ) : u { € C 2 (T) 


u 


^ k (J * q k ,s) + qp,s 


* u 


k=—d 


k=—d 


and 


u\e)&0 = 0, 



(2.13) 


(2.14) 


Due to the invariance by rotation of the model (11.71) . is linked to Lq s in an obvious 
way: L^^u^(-) = L 0i 5 tt(.), where u^(-) = u( ■ — ip), so that the operators {L^^)ipG T 
obviously share the same spectral properties. For any operator L, the usual notations 
a (L) (resp. p(L) and R(X,L)) will be used for the spectrum of L (resp. its resolvent set 
and its resolvent operator for A € p{L)). 
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One can prove (see [21], Theorem 2.5 and Appendix [Bl below) that for all 0 ^ <5 ^ SfiK), 
L^g is closable in Hff 1 , sectorial, has 0 for eigenvalue, associated to the eigenvector dgq^^g, 
which belongs to the tangent space of M in q^ g (this reflects the fact that the dynamics 
induced by (II.711 on M is trivial) and that the rest of the spectrum is negative, separated 
from the eigenvalue 0 by a spectral gap 7 l > 0. More details about these questions are 
given in Appendix iDl 

The fact that the eigenvalue 0 is isolated from the rest of the spectrum a(L^^g) \ 
{0} implies that Hfi 1 can be decomposed into a direct sum T^ t g © N^ t g, where T t ^g = 
Span^dgq^j) such that the spectrum of the restriction of L,^g to N^ t g (resp. T^ t g) is 
\ {0} (resp. {0}). We denote by P®, s the projection on along and 
g = 1 — P®g- Both P® s and P^ g commute with L^g- In particular, for all ip E T, 
5 > 0, there exists a linear form satisfying, for all u € Hfi 1 

P °,6 U = Vi>,s(u)d e q^,g . (2.15) 


We also denote by Cp and Cl positive constants such that for all u € H d \ t > 0 : 

\\ P $,5 u \\-l,d ^ Cp\\u\\-i t d , 

\\ P $,,6 u \\-l,d < Cp\\u\\-i, d , 


\e tL ^Plgu\\_ ld < C L e-^\\Plgu\\_ ld 


\e tL *' s u\ 


-1 4 


< C L 1 + 


Vi 


\u\ 


-2 ,d ' 


(2.16) 

(2.17) 

(2.18) 

(2.19) 


Inequality (|2.18l) is a consequence of )24| . Theorem 1.5.3, p. 30 and (12.191) is proved in 
Proposition IB.71 in Appendix [Bj Once again, we will often drop the dependency in the 
parameters ip or S in P^ s and P | s for simplicity of notations. 

A consequence of the contraction (12.181) along the space g is that M is locally stable 
with respect to the evolution given by m (see for example exercise 6 * of the Chapter 6 
of (24) . or Theorem 2.2 of [21] for our particular model): for any po i n a neighborhood of 
M, there exists ip € T such that the solution of (11.71) converges to q^j exponentially fast 
(with rate given by 7 l). 


2.4. Dynamics of the empirical measure. The starting point of the proof of Theo- 
rem l2.3l is to write the semi-martingale decomposition (see Proposition l3.1D of the difference 
between the empirical measure defined in (11.61) and any element of q^j E M. Namely, 
define the process 1 1 ->- vgr t , t ^ 0 by 

v N,t := A l N,t ~ V,Si * = ~d, ■ ■ ■ j d. (2.20) 

The point is to write a mild formulation of this semi-martingale decomposition that makes 
sense in the space Hf 1 (recall that /xjv,t and vN,t belong to Hfi 1 due to (12.41) 1. This mild 
formulation involves in particular the semi-group e tL, t’ s of the operator L^g (12.131) so that 
one can take advantage of the contraction properties of this semi-group in the neighborhood 
of the manifold M. 


Proposition 2.7. For all K > 1, for all 0 ^ 5 ^ d(K), the process (i7v,t)t> 0 defined 
by (12.201) satisfies the following stochastic partial differential equation in C([0, + 00 ), Hf 1 ), 
written in a mild form: 

VN,t = e^^z/jv.o + [ ed -s ) L ’/ , .' s ( D^ — dgRiy(vN,s)) ds + ^N,t, N ^ 1, tfi 0, (2.21) 

Jo 
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where 


Dn = D N tf >( 5 := -de ( q^,s ^ (A k N - A k )(J * q^ s ) 


k=—d 


( 2 . 22 ) 


\k=—d 


Rn(vn,s) = Rn,iP,5{^N, S ) ■ = ( ^2 (X k N - X k )J *q§A u N,s+q^,5 (^N-^ k )i J *^N,s) 

J k=-d 

+ ( '^2 A k N J * \ un , s , (2.23) 

\k=-d J 


and Z/v,t is the limit in H d 1 as t' /*• t of ZN,t,t' defined by 


z N , t ,t'(h) = ^2 w J2j ' 

i=—d j =1 J U 


e^+h 


(ri(s))dB)(s), (2.24) 


that we denote 


d \i ft 

= E jvtE/ 8 > 

i=—d j=l Ju 


f d 0 


Jo 

\ J 


(^.(s))dS;( S ), (2.25) 


and where all the terms in (I2.2ip make sense as elements ofC([0, oo), H d ). 


The proof of Proposition 12.71 may be found in Section [31 The term Zjy,t in ( 12 . 211 ) 
represents the effect of the thermal noise on the system. The term involving Dn is the 
one that produces the drift we are after on the time scale IV 1 / 2 ^ when the empirical 
measure /ijv.t is close to the manifold M. To make this drift appear, we rely on an 
iterative procedure, as explained in Section [2T6l 


2.5. Moving closer to the manifold M. We place ourselves in the framework of The¬ 
orem [232 we fix Co > 0 and suppose the existence of a probability measure po € Hf 1 such 

that dist^-i (po, M) ^ eo with P ^||^at,o — Po||_i^ ^ 0 as IV —>• oo, for all e > 0. 

The constant eo will be chosen small enough in Section El 

The first step in proving our result is to show that the empirical measure p,N,t reaches a 
neighborhood of size iW 1 / 2 in a time of order log N. We use the projection defined in the 
following lemma, whose proof can be found in Appendix [Cl along with several regularity 
results. 

Lemma 2.8. There exists a > 0 such that for all h such that dist^-i (h,M) ^ a, there 

exists a unique phase if =: proj^(/r) € T such that P®(h — q^) = 0 and the mapping 
h i-a proj M {h) is C°°. 

From now on, we fix a sufficiently small constant (, more precisely satisfying 

C < l • (2.26) 


We prove the following result: 
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Proposition 2.9. Under the above hypotheses, there exists a phase 9 q € T, an event B N 
such that P (B N ) —> 1 and a constant c > 0 such that for all £ > 0, for N sufficient large, 
on the event B N , the projection i/’o = f-’o = proj M (^Af.w 1 / 2 ^) well-defined and 


Un,n 1 / 2 t$ 


Qi/jo 


^ AT-V2+2C 

-1 ,d 


(2.27) 


and 

\ifo-0 o \^e, (2.28) 

where t(f = cfV -1 / 2 log IV. 


We refer to Section [6] for a proof of this result. Since it relies on a discretization scheme 
similar to the one we introduce in the next paragraph, we leave the details to Section [6] 


2.6. Dynamics on the manifold M. We now place ourselves on the event B N (see 
Proposition 12.91) . so that on the time N 1 //2 fj^ we have WfJ-N NUHff ~ Ift/’oll-M ^ iV -1 / 2+2 ^ 
where ifo = proj M (/% n 1 / 2 ^) • The point is to analyse the dynamics of (12.2111 on a time 

scale of order IV 1 / 2 , using the knowledge we have on stability of the manifold M (recall 
(11.121) ). The following iterative scheme we introduce is similar to ones used in |7, .12] . 


The iterative scheme. We divide the evolution of the dynamics (12.2111 in time intervals 
[T n ,T n+ i] with T n = N 1 / 2 ^ + nT where T is a constant independent of N, satisfying 
T Js 1 and 


e" 7iT sC 


1 

4 C L C P ’ 


(2.29) 


where the constants Cl and Cp where introduced in Section T2. 31 The number of steps np 
is chosen as 


n f 



(2.30) 


The intuition of this discretization is the following: if for a certain n = 0,1,..., np — 1, 
the process p,p n = hN,T n is close enough to the manifold M, we can define the phase a n 
of its projection on M by: 

a n ■= Woj M (UN,T n ) ■ (2.31) 

This projection is in particular well defined when \\/J-T n —Qa n -i ll-M ^ where the constant 
a > 0 is given by Lemma 12.81 

To ensure that the process does not escape too far from M , we introduce the following 
stopping couple (where the infimum corresponds to the lexicographic order): 


(n T ,r) = inf{(n,f) € {1,... ,np} X [0,T] : ||/iT„_i+t - q an -i \\-i,d > cr} . (2.32) 

Using ()2.32f) . we can define the following sequence of stopping times (r n ,n = 1 ... np): 


( T if n < n T , 
| r if n ^ n r , 


(2.33) 


and consider the stopped process li( n An T -i)T+tAr n - The projection of this stopped process 
is well defined on the whole interval [Tq. T nf \, so that we can now define rigorously the 
random phases i/’n-i defined as 


ifn— 1 ■ P^®j M (U(nAn T — l)T+tAr n ) > ^ l,...,np. 


(2.34) 
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ip n -i corresponds to the phase of the projection of the process ji unless it has been stopped 
and in that case, it is the phase at the stopping time. The object of interest here is the 
process u n ^ defined for n = 1 ,..., np as 

Vn,t ■— M(nAn T — l)T+t/\r n Qipn-l • (2.35) 

Using (12.211) . we see that this process satisfies the mild equation 

g (iAr s)L^ n _ 1 _|_ R x j Jn _ 1 (y n>s ')') ds + Z n ^t/\ T n , (2.36) 

where D^, nl := D N ^ n _ uS (recall (12.221) ). R^flvn) '■= Rn^^^u) (recall (I2.23D ) and 
Z n j is defined as 


'n,t 




0 


rtA 

Jo 


ZnM =Y.NlY.i 8 » 

i=—d j =1 


/0 


e ( * h 


(<Pj(T n -i + s)) dB*-(T n _! + s). (2.37) 

Note that we drop here the dependence in N and 6 for simplicity. 


Controlling the noise and a priori bound on the fluctuation process. A key point in the 
analysis of (|2.36l) is to show that one can control the behavior of the noise part Z n t in 
(I2.36p along the discretization introduced in the last paragraph. More precisely, for £ 
chosen according to (|2.26l) and some positive constant Cz and defining the event 


A n = A n (C z ) 


sup sup \\Z n fl 

l^n^rifO^t^T 


1 ,d ^ Cz 



(2.38) 


the purpose of Section 0] is precisely to prove that P(A Ar ) tends to 1 as N —>• oo. With the 
knowledge of (12.381) . one can prove that the process v n remains a priori bounded: using 
that the sequence of the disorder ^ i is admissible (recall Definition 12. ip . we prove in 
Proposition 15.11 Section [5l that on the event A N n B N , 

sup sup \\v n ,t\\-id = 0(N~ 1/2+2 () , (2.39) 

1 < n ^ rif te[0,T] 


as N —> oo. 


Expansion of the dynamics on the manifold M. The last step of the proof consists in 
looking at the rescaled dynamics of the phase of the projection of the empirical measure 
on M, that is the process 

E t N := VW , (2.40) 

where (ip n )o ^n^n f is given by (|2.34|) and 


n t ■■= 



(2.41) 


Namely, we prove in Propositions 15.21 and 15.31 that, with high probability as N —» oo, the 
following expansion holds: 

E t N = if 0 + b(£ N )t + Q(N -V 4 + 2 C), (2.42) 


where b is the linear form of Theorem 12.31 and that ji N N 1 / 2 t i s close to q^N with high 
probability. 
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2.7. Organization of the rest of the paper. Section 0 is devoted to prove the mild 
formulation described in Paragraph s. 41 The control of the noise term in (12.371) is addressed 
in Section [4J The dynamics on the manifold M and the approach to the manifold are 
studied in Sections [5] and [6] respectively. The asymptotics of the drift as 5 —>• 0 is studied 
in Section 17.21 We compile in the appendix several spectral estimates and expansions in 
small 5 used throughout the paper. 


3. Proof of the mild formulation 


Define Lq d as the closure of the space of regular test functions (u d ,..., u d ) such that 
f T u k = 0 for all k = —d,..., d under the norm 


o ,d 



(3.1) 


and the space H d (a ^ 0) closure of the same set of test functions under the norm 
(denoting || • ||q the L 2 -norm on T) 


1/2 


I a,d ' 


(~A,) a/2 i 


0,d 


E (~ A ) 


a/2 u k 


\k=—d 


(3.2) 


where denotes the Laplacian on T 2d+1 . We denote by H d a the dual space of H d . We 
also write, for any bounded signed measure m on T, the usual distribution bracket as 

(m, /) := f f( 9 )m(d 0 ), 

Jt 

and for any vector (to 1 ,..., m d ) of such measures 

m,F) d := Y j \ i {m\F i )= ^ V 


F i (6)m i (d9), 


the corresponding bracket weighted w.r.t. the disorder. Obviously, when the above mea¬ 
sure coincide with an L 2 function, this expression coincides with the L 2 scalar product 
(•, -) 2d associated to m- 

This section is devoted to prove Proposition [221 We begin first with a weak formulation 
of the SPDE (1T2TD . 


Proposition 3.1. For all K > 1, for all 0 ^ 6 ^ <5(li) ; f or an V (t, 0) FfiO) = 
(F t ~ d (9),...,F t d (0)) eC 1 ’ 2 ([0,+oo) x T,M) such that f T F t (9)dd = 0, 

(vn,s ; d s F s + L^ d F s ^ d ds + J (Djy , F s ) d ds 

+ [ (Fn^n^) i dgF s ) d ds + Tf y i , N ^ 1, t ^ 0 , (3-3) 

J o 

where Djy, Rn(vn) are respectively defined in (12.221) and (|2.23l) and 

d N* 

M N,t ■■= 

i=—d j=l 


S' 

Jo 


d e Fi{Ms))dB){s) 


(vN,t,Ft)d = (vn,OiFo 


)d+ f 

Jo 


(3.4) 
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In (13.31) . the operator L*^ 5 is the dual in of the operator in (12.131) : 


1 d 

(L^y = -^dgv 1 + dufdev 1 + (dev 1 ) ^ X k J * qfy s 


k=—d 


(dev*) J 2 A k J*qis) dd 

\ k=~d / 

d 

k J*(4,sdev k ). (3.5) 


k=—d 


We refer to Appendix [B] (see in particular Propositions IB.31 and IB.41) for a detailed 
analysis of the spectral properties of the operator and its dual L ^ 5 . All we need to 
retain here is that when 5 is small, the operator L^g is sectorial in HI 1 and generates 
a Co-senri-group t H > e tL ^’ s in this space. Moreover, on the space Lq d , one has that 

( e tL ->b s)* = e tL *’ 5 . Since the phase -0 is not relevant in this paragraph, we write for 
simplicity qg, L$ instead of q^ t g and L^ j( y. 


Proof of Provosition \3.1\ Note that, using the definition of </(•) in (11.41) and of the empir¬ 
ical measure HN,t hr (H-61) . the system (11.51) may be rewritten as 

d 

d tpj(t) = 5u> 1 dt + ^ \ k N J * n k (<Pj(t)) dt + dBj(t), i = —d,...,d. (3.6) 

k=—d 

Consider (t,6) F t (6) = (Ft(^))i=i,—,d € C 1,2 ([0,+oo) x T, M) rf such that for all t ^ 0, 

f T F t (9) d 9 = 0. An application of Ito Formula to (11.51) gives, for i = 1,..., d, j = 1,..., N l , 
t> 0, 


FI 


M(t)) = F*(<p)( 0 ))+ fdsF&iWds + l f d 2 e Fl(v)(s))ds 

JO 1 Jo 

+ doFs(p l j(s)) (du l + ^ X% J * ( s ))\ ds + deF):(ip l j(s)) dBj(s) 


k=—d 


After summation over j = 1 ,,N l , we obtain, for i = 1,... ,d, 
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Replacing fi l Nt by v l Nt + q\ in (13.71) (recall (12.201) 1. we obtain 

{ v N,ti Ft) + (Qh Ft) = { u n,qi K) + (vb Fq) 

+ J (^N,s + QS’^sFl + -dgFg + d eK + Y^ A k N J * {v% )S + ^ ds 

JV* ^ 

+ 4E / 9»A>5(»))dB‘(«) 

3 =1 ° 

= {F nq , Fq) + (\i' l iy ! si &sFl + - dgFg + d e Fl + Y^ A k N J * q k ^j ^ ds 

+ ( ,y k T ,s- doFg Y, A)y («/ * u N,s) \ ds + (qb dgF ,l s Y^ A^( J * v n,s) \ ds 

+ (qi Fq) + j* (qi d s Fi + l -d 2 e Fl + d e Fl ^ + ]T \ k N J * q^j ^ ds 

jV* i 

+ 4E/ 6eFt(v)(s))dB‘(s). (3.8) 

3=1 J ° 


Since by definition q$ is a stationary solution to (|1.7D . one easily sees that 


\,F{) = (qb Fq) + / (qld s Fi)ds, i = l,...,d, 0 ; 

Jo 


(3.9) 


and 


0 = l^-dU - Su'doqi - dg ^ fl Xkj * > F^j 

= lq\} T dlFl + 5Fd 0 Fl + d 0 Fl Y, A k J*q%) • 


(3.10) 


k=—d 


Summing (13.8|) over i = —d ,..., d and using (13.9p and (13.101) . we obtain 


(z 'N,t, F t )d = Wn, o, F 0 ) d + / / v N>s , d s F s + -dgF s + 5d e F s <g> w + d e F s Y x n j * Qs ) ds 


k=—d 


k=—d 


+ j \VN,s,9gF s Yj Ajv</ * fN,s J ds + J (^qs,dgF s Y^ ^n^* u n,sJ ds 

+ T (qsY^- Xk ) J *^^Fs\ d s + M^ t , (3.11) 


k=—d 
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where M^ t is defined in (13.41) and where we have used the notation F <S> u = 

Note that 

/ d \ d d 

U,d e F s ^2x k J^ k N>s \ VA k (qi,deF' s J*v k NiS ) 

\ k= 1 / d i=l k=—d 

d d 

- - V V \ i \ k (v' K _„J,(q}deFl)) 

:—d k=—d 

v s ,f^\ k J*(q k d e F k )\ . (3.12) 

fc=-d / d 

The result of Proposition 13.11 is a simple reformulation of (13.111) using (13.121) and the 
definition of L d in (13.51) . □ 

We are now in position to prove Proposition 12.71 

Proof of Proposition \ 2. 7\ Let us apply the identity (13.31) of Proposition 13.11 in the case of 
test functions F s of the form 

F s = e (t - s)L sh, 

for any test functions h of class C 2 on T. Then d s F s = —L d F s and one obtains 

C VN,t,h) d = (vN,o,e tL * 6 h) d + j (p N , ds + J (Rn(vn,s) , d e e^~ s)L sh^ ^ ds 

+ • (3.13) 

We aim at proving that one can write a mild version of this weak equation and that this 
mild formulation makes sense in Ftf 1 - Consider a sequence (vi)i^i of elements of Lg d 
converging as l —>• oo in H d 1 to vn.o € Hff 1 . Then, for h of class C 2 , 

(vi, etL * sh ) d = ( v h etL * Sh ) 2d = ( etLSv F h ) 2 , d = ( etLSv i ) h ) d • (3.14) 

By continuity of e tLs on H d l , e tLs vi converges in Hf 1 to e tLs VN, 0 i as l —>• oo. In particular, 
for all h € H d , 

\(e tLs VN,o , h) d - {e tLs vi , h) d | SC ||h|| M \\e tLs v N ,o ~ e tLs vi\\_ 14 ~>i-* oo 0, (3.15) 

so that, at the limit for l —>• oo, for all t ^ 0, 

(v N , o, etL * sh ) d = ( etLs "N, o, h) d . (3.16) 

Since the function Dn defined in (12.221) is regular, it is straightforward to prove in the 
same way that 

(d n , e^~ s)L ^h}j i = (e^b'Djf , h) . (3.17) 

The continuity of the mapping t e- e tLs VN,o and t e^~ s ^ Ls Dn ds in H d l is immediate 

from the continuity of the semigroup in HJ 1 . 

We now focus on the term Rn{vn,s) defined in (|2.23l) : consider ('w s ,i)i ^ i a sequence 
of elements of Lg d converging in H d 1 to vn,s (consider for example w s j = 4>i * vn,s for a 
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regular approximation of identity ((j>i)i ^ 1) and define 


\k=—d 


k=—d 


Rs,l ■= ( A V * 9* ) + QS iz ( X N - A ")( J * A,s) 

+ ( ^ A^rJ * A,s\ w s,i ■ (3.18) 

\k=-d ) 


2 A 


(3.19) 


For any l ^ 1, the following identity holds: 

(r s ,i , doe( t ~ s ' >L s hj^ = (^R s ,i , = - (e^ 8 ^ 5 d g R s ,i , /i^> 

Since h is regular and R s j converges in H d l to Rn{vn,s), the lefthand part of the previous 
identity converges as l —>• oo to {Rn(^n,s) , d g e^~ s ^ L s h) Moreover, for all h regular, using 
the estimate (IB.24|) on the regularity of the semigroup e tLs , (note in particular that 
can be extended to a continuous operator to HT 2 to LT7 1 , see Proposition IB.71 below). 


e (t s ^d g (R St i - R n (v n ,s)) , h 


€ 


^ C 


< c 


1 ,d 


e (t s '> Ls d g (R St i — Rn{vn,s)) 

1 


1 + 


—l,d 

\\d g (R s ,i — Rn(vn, 


,d ’ 


y/t — S 

, d “ R n(vn,s) l|_ M • 


Since the last estimate is true for all h regular, one obtains that 

1 


e h s ) L ^d g (R s ^ - Rn{v n , s )) 


-l a 


^ C 1 + 


y/t — S 


II Rs,l ~ RN(vN,s)\\_id- (3-20) 


Since R s j converges to Rn(^n,s) in H d 1 , one obtains that one can make l oo in (13.191) : 

Rn(vn,s ), d g e (t ~ s)L sh)^ = - (e (t ~ s)Ls d g R N (v N ,s ), hj ^ . 

The same argument as before shows also that 


e h s ' ,Ls d g RN(vN,s) ^cfl+ , ||-Kiv(kv,s)IL 1 d 

-1 ,d \ y/t — S J ’ 


< C 1 + 


y/t — S 




sC Ctt 


\ 


(3.2D 


k=—d 


where we used (|2.4I) . The inequality (13.211) implies that the integral Jg ||e^ s ^ Ls dgR^i^N^W^ 
is almost surely finite. Using @3], Theorem 1, p. 133, we deduce that fg e^~ s ^ Ls d g RNiyN,s) ds 
makes sense as a Bochner integral in HJ 1 . The continuity of t *->• fg e^~ s ^ Ls dgR^^N,/) ds 
in H/ 1 is a direct consequence of the bounds found in Proposition IB.7I 

It remains to treat the noise term in (|3.13l) . The precise control of this term is made in 
Section [4] below (see in particular Proposition 14.11) . We prove actually more in Section [4] 
since we have to take into account the dependence in N, which is not important for 
this proof. Let us admit for the moment that the proof of Proposition 14.11 is valid. In 
particular, one deduces from (14.281) and an application of the Kolmogorov Lemma that 
the almost-sure limit when t' /* t of t # defined in (12.241) exists in H 7 1 . The continuity 
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of the limiting process t i-a Z^,t in H, 1 comes from (14.3p and the rest of the proof of 
ProDosition l4.ll It is immediate to see from (13.41) that for all F regular, (Zjv jt , F), = 
where we see the term Z n> t as a vector (Z~f,...,Z!* t ), with Z^ t {h) = pr Z ni t(hk) and 
hk = (0,..., 0, h k , 0,..., 0). One concludes from everything that we have done that, for 
all h regular that 

{vN,t, h) d = (e tLs u Nfi , h) d + (e [t ~ s)Ls D N - e ^~ s)Ls dgR N (i/ NtS )') ds , hj 

+ (Z Ntt ,h) d , (3.22) 

where everything above makes sense as element of HT 1 . Since this is true for all h regular, 
the identity (12.21[) follows. Proposition 12.71 is proved. □ 

4. Controlling the noise 

This section is devoted to control the noise term Z n>t defined in (12.371) . More precisely, 
we prove the following proposition (recall the definition of A N = A N (Cz) given in (12.381) 1. 

Proposition 4.1. For all £ > 0, there exists a constant Cz such that P(A^) -A 1, as 
N -A oo. 

To prove Proposition 14.11 we rely on the two following lemmas: 

Lemma 4.2 (Garsia-Rademich-Rumsey). Let x an d 4' be continuous, strictly increasing 
functions on (0, oo) such that %(0) = 'P(O) = 0 and liirp'l'(t) = oo. Given T > 0 and 
(f continuous on (0, T) and taking its values in a Banach space ( E , ||.||), if 




to Jo 


-*0011 


X(l t ~ s|) 


dsdt ^ B < oo, 


then for 0 ^ s ^ t ^ T, 


~ *0011 ^8 f T 1 

Jo 


) x(du) • 

u z 


(4.1) 


(4.2) 


, Theorem 2.1.3. The second result estimates 


Proof of Lemma 14.21 mav be found in 
the moments of the process Z nt \ 

Lemma 4.3. For all e > 0 and all integer m > 0, there exists a positive constant C m> 
such that for all 0 ^ s < t ^ T, 

c 


— Z r 


|| 2 m 

-,s II —l,rf 




N m 


'-({t- S ) m ( 1 /2-2 £ ) + _ 

U \ / 


(4.3) 


Let us first prove Proposition 14.11 relying on these two lemmas. 

Proof of Proposition \f~l\ Using Lemma 14.31 we can apply Lemma 14.21 with the choices 

(fit) = Z n!t , xi u ) = and T(ii) = u 2m , (4.4) 

which implies that there exist a constant C (depending in m, £ and £) and a positive 
random variable B such that for every 0 ^ s < t ^ T: 


Zt! I Zr,, 


12m. 

1-1, d. 


< C{t-sfB, 


where B satisfies 

E (B) < 


C 


ff 


t — s | m ( 1 /2-2e)-2-C _)_ u _ s im—2—C 


ds dt . 


(4.5) 

(4.6) 
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A simple integration shows that E(£>) ^ (T m ( 1 / 2 - 2£ )-( + T m ~^), whenever m{ 1/2 — 

2s) — ( > 1 and m — ( > 1, that is when m > We can fix for example s = 1/8 and 

choose an integer m such that m > 4(1 + £). Since T ^ 1, we have E (B) ^ and 

we obtain: 


7 _ y II2 m N 

^n,t ^n,s\\-id 


E sup ,, 

\0^s<t^T \t ~ Sp 


^ C- 




N r 


(4.7) 


which implies 


P sup 

\ o ^ t ^ T 


We deduce 


\Z„,,U4 > \/J WC ) < ( o sup t l|Z„,,ll5,u) 

( 


N m 

< --AT- 2m CE 

^m—Q 


sup ■ 

V 0 ^ t < T 


y 11 2m \ 

] < CN~ 2m<: . (4.8) 


sup sup 

l^n^rifO^t^T 


\Z n ,t\\-i, d > < Cn f N~ 2m <, 


(4.9) 


which tends to 0 as N —>• oo if we choose m > since n/ = 0(N 1 / 2 ). Proposition 14.11 is 
proved. □ 

Let us now prove Lemma 14.31 

Proof of Lemma E3 Our aim here is to get the appropriate bounds for the process Z. 
We follow mostly the ideas of [21] . Recall that 


d \i Ni rTn-i+t 


Z n ,t{h) - N i 


de 


i=—d j=l jT n-l 


(■ t-u)L* , \ * 

e h 


{^{u))dB){u). (4.10) 


Let us define the process Z n ^t' for 0 <t'<t as 


d fffl_ rTn-i+? 


z UAh) = £ wT. 


N' 


dfi 


i=—d j = l J Tn — l 


M u > L in-i h y (pi(u))dBj(u). (4.11) 


Our aim is to estimate for 0 < s' < s < t, s' < t 1 < t and for all integers m > 0 the 
moments E(|| Z n ^ — Z nss i ||^ d ). We can decompose Z n t t i — Z ns s i as follows: 



^ n,t,t ' 

where 


M n,s',s,t( h ) 

d x i Ni r Tn - 

= T,w-Z,J t , 

i=-d j'=l n - 1 

and 



Zn,t,t' Z n s s t — M n s i s t + M n)S ij.ij. , 


de 


L \i Hi rTn- 1 +t' 


(t—u)L*, (s—u)L*i 

’ Vn-l — g V ' Vn-1 


(4.12) 

(¥>}(«)) d5*(u), 

(4.13) 


= E 


*-* N- ‘-- /T-I l J 

i=-d j = l^^n-l+S 


cb 


e ( * (^.(u))dR;.(u). (4.14) 
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and using Holder and Burkholder-Davis-Gundy inequalities, we obtain for the terms in¬ 
volving M 1 


E lJ2\<s',sA h l)f 


1=1 


^1 5^2 — 1 

oo 

OO 

^ C m ^ E [M^ )Si ,(hq)] s/ ... E K. !S ^,J] s/ 

^1 5/2 v*5^m = l 




— I d N1 rTn-i+s' . . 2 

g. E EE/ t 


h,l2,...,l m =l \i=—dj=l jr ^ n 1 


E dir 


jVi /•T„_ 1+S ' 2 

E(<"#/j) dir 


££ Jr 

— 1 J J- n —1 


i=-dj=l 

00 d N i <.T„_!+s 


c ™ eee/ t 

\Z=1 i=-dj=1 J 1 

( d Ni rT n -i 

cJee/ t 

\i=—dj =1 J lrL ~ 1 


e ( 


dtt 


E||<tf M || 2 _ M d« 


(4.21) 


The same work can be done for the terms involving M 2 , which leads to 

d ^ /•T n _ 1 +s' 


12m 
^n,s,s' II—1 


E Z„ t t > — Z, 


E||<"rll 2 -Mdu 

d ^ r T n-l+t’ 


l d Ni rT n -i 

cq E E J T 

\i=—dj =1 J dn ~ 1 

( d L fJn-ltr 

EE/ E||H 2 "e M dn| . (4.22) 

i=-d 7 = 1 


It remains now to hnd appropriate bounds for E and E||C/;^|| 2 _ m . On 

hand, for h E H\ . since (5m € H~ 1 l 2 ~ e for all e > 0, we have 


r 2 4d II2 


one 




_c 

To 


& 


Xt-^h 




l/2+e,d 


c_ 

N 1 




3/2 +e,d 


(4.23) 


For the rest of the proof, we set e = | (any e € (0,1/4) would be sufficient). Applying 
Proposition IB.61 with (3 = 1/4 + e/2, we obtain, for any 0 < 7 < 7^*, 

lO*)l 4( W - «)-*'«-/*) ||fc|| w , (4.24) 

which means that |l/ 2 ’/ J t ||-i, f i ^ (l + e _7 d _ “)(t _ u) -1 / 4-6 / 2 ). On the other hand, 
proceeding as before, we get the bound: 




c 


N l 


( e (t " s)L * - 1 e (s -“ )L ^y 


3/2+e,d 


(4.25) 
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Applying Proposition IB.6I with f3' = 1/4 + e/2 and (3 = 1/4 — e, we get for all h € IL 2 £ , 
\ l[e (ts)L } _i]^|| 3/2+M ^ C £ (t — s) 1 / 4_£ ||(1 — . (4.26) 

For 5 = e^ s ~ u ^ L sh and using again Proposition IB.61 with this time /3 = 1/2 — e/2, this 
leads to 

(4.27) 

which means that ^ — s) 1 / 4_e (s — u)~ l / 2+£ / 2 e ~ 1 ^ s ~ u " > . We can now 

estimate (14.2211 : using that N ^ ciV* ^ CN , we obtain 



r^m 

< s ) m(1/2 ~ 2£) + (t' - s') m(v2_£ ) + (i' - s') m ) • (4.28) 

Taking t' /* t and s' s and using Fatou Lemma, we deduce the result. □ 

5. Dynamics on the manifold M 

The purpose of this section is to prove the results described in Section 12.61 concerning 
the process v n defined in (|2.35l) . 

Recall that the scheme defined in Section [2761 starts at a time t\f = 0(N~ ll/2 log N), 
such that there exists an event B N with P (B N ) -A 1 such that on B N if we denote 
ipo = proj jvf (/r jy jyi/2^ ) then \\n NN i/ 2 t N - q^ 0 \\-i,d ^ N~ l t 2+2C -. In other words, the 
initial condition of the scheme satisfies ||z^i,oII—i,rf ^ N~ 1 ^ 2+2 ^ on R^.The existence of 
these times and event B N will be proved in the Section [6j The first result proves 
estimate (12.391) : 

Proposition 5.1. There exists an event with P(D) V ) -A 1 as N -A oo such that, 
almost surely on , 

SUp SUp ||z^rx,i||-l,rf = 0(iV _1/2+2C ), (5.1) 

1 ^ n ^ nf te[0,T] 

where the error 0(N~ 1 ^ 2+2< ^) is uniform on fb/ . 

Proof of Proposition 15.11 Recall the definition of the event A N in (12.3811 and define ii] 2 := 
A n fl B N . Since the purpose of Section |4] was precisely to prove that P(A iV ) -A 1, we 
obviously have that P(fl) v ) -A 1, as N -A oo. 

Throughout this proof we work on the event and proceed by induction. We already 
know that ||ui,o||-i,d ^ _/v _1 / 2+2< ’. If we suppose that ||iAi,o||-M ^ _/v~ 1,/2+2 ^, then from 
the mild formulation (I2.36p . from (12.1811 and (12.1911 and from the estimates on the noise 
term Z n)t on C A N , we obtain 

IK«||-M < C Le ~^" 

+ 2TCl\\D^ ti _ 1 II _l,rf + Cl(T + 2T 1 / 2 ) sup \\R^ n _ 1 {l / n,s)\\-l,d 

0 ^ s sC t 

+ T l/2 N -l/2+( ' ( 5 . 2 ) 
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Since the sequence (uii)i ^ i is admissible (recall Definition 12.11) . we have 

< CN- 1 / 2 max |^| < CiT 1 / 2 ^. (5.3) 

k=—d,...,d 

Define the time t* as 

t* ■■= inf [t € [0, T] : |K, t ||_ M > 2C L N~ l / 2+2 ^ } . (5.4) 

Obviously t* > 0 and if t ^ t*, one readily sees from (12.2311 that 
SU P 11 #^—i(Di,s)||-m 

0 ^ s ^ t 

^ cf sup \Wn,s\\-i d + N ~ 1/2 max |£&| sup |K,s||-m^) < CiV _1+4C . (5.5) 

Vo^s^t k=-d,...,d o^s^t ) 

Putting together (15.211 . (15.31) and ()5.51) gives that t* = T if N is large enough. Consequently, 
by construction of the stopping time r n in (12.33[) . one has that r n = T and the choice of 
T (recall (12.29[) 1 implies that 

\Wn,T\\-i,d < t-^- lV _1/2+2? . (5.6) 

ZO p 

To conclude the recursion it remains to show that HiAi+roII-M ^ -/V _1 / 2 + 2 t\ To do this, 
let us write v n +\fi in terms of 

Di+i,o = Qiji 7 !_i + v n,T — Qil> n ■ (5-7) 

Since P^u n+ i,o = zAi+i,o : where we recall that P^ is the projection on the space N^ n , we 
can rewrite it as 


Dl+1,0 — PlpniQlpn-l T Vfl,T Qlpn) 

= P$ n (.Qil>n-l - qi/in) + (P$ n - + P^> rl _ 1 Di,T ■ (5.8) 

Since q^ n _ x — q^ n = (V’n-i — V’n)^ + 0((ip n — ip n -i) 2 ) (and this estimate makes sense 
in H J 1 ) and P^ n dgq^ n = 0, the first term of the second line of ()5.81) is of order 0((ip n — 
ip n - 1) 2 )- Using the smoothness of the projection proj M (Lemma 12.811 . 

IV’n — V’n-ll = I P r °jAf (AhnAn T -l)T+iAT«) ~ P ro jiW(^((n-ljAnT-ljT+tAr"- 1 )l 

U C*11A^(ri.Arj.T- — 1)'Z~'+iA t - ” A^((n— l)An T — l)T+t/\r n ~ 1 II — l,d 

^ C\\v n - hT \\-i,d + C\\v n - lfi \U >d < CUV-V^. (5.9) 

Combining the last two arguments, we obtain that the first term of the second line of (|5.8I) 
is of order 0(iV~ 1+ U). For the second term, the smoothness of the mapping ip i-a P*, gives 


||(P^-P^_> n , T ||_ M < C\lP n -lPn-l\\WA-l4 < CN~ 1+ <. 
Taking the H^ 1 norm on the two sides in (15.8}) . we obtain 

IIPn+1,01| — 1 ,rf ^ l|-f > ^„_ 1 Z'ra,T||- 

which implies the result for N large enough. 


i, rf + 0(iV" 1+4C ) < IiV- 1/2 + 2 C + O(7V- 1+4 0, 


(5.10) 

(5.11) 
□ 


We are interested in the rescaled dynamics of the phase of the projection of the empirical 
measure on M and in particular use the rescaled discretization of this phase dynamics given 
by the process (recall (12.4211 i. 
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Proposition 5.2. There exist a linear form b : M 2 ^ 1 —>• M and an event 12^ satisfying 
P (£2^) —>• 1 as N —>• oo such that on the event £2^ we have for t £ [to ,tf]: 

*t N = V>o + b(£ N )t + 0(N - 1 / 4+2C ), (5.12) 

where the 0(A r ^ 1 / 4+2 ^) is uniform on f2^. 

Proof of Proposition 1 5.121 We work for the moment on the event 12^ defined in the proof 
of Proposition 15.11 Using Proposition 15.11 Lemma 1C.II below and the fact that if n = 
proj+ u n,T), we have the following first order expansion of P f N in (15.1211 (recall 
the definition of p in (12.1511 and note that there are 0{N 1 ^ 2 ) terms in the sum): 


:= V>o + Ep b„- 1 K,T) + o(iv- 1/2+4C )- 


(5.13) 


n=l 


Let us now decompose the term p^„_ 1 (^n,T) ) using the mild formulation (12.361) . Re¬ 
mark that p^, n _! (e 1 **- 1 u n fi) = p^i„_i (zAi,o) = 0 and that p^„_ x (e (t_s ' |L ^- 1 = 

Pp n _i(A/> n _i)- Note that Proposition 15.11 shows that r n f = T on 12^, so that the time 
integration in the mild formulation (I2.36P does not involve any stopping time. Hence it 
remains, since D^ n _ 1 has no dependency in time, 

Pp„_iK,r) = %_i(%_j) - J p^ ^ e (t -* )L *n-i doR^^iy^sfj ds + Py,^ (Z n , T ). 

(5.14) 

Using ()2.19f) and (15.511 


PV- 


_1 ^ e (t S ) L V>n-l cl 0 i?,,/, ri _ 1 (v n , s )) ds 



r T 

j ds 

< 

Jo 


^ c 


f 

Jo 


e (t s ^ L,l ’ n - 1 dgR^ n _ 1 (un tS ) 

1 


—l,d 


ds 


1 + 


< C(T + Vf)N ~ 1+4C , 


II R il>n-A u n,s)^_ l d ds 

(5.15) 


which leads to 


PP„_:lK,t) = +p^_i(2„,t) + 0{N 1+4C ). (5.16) 

We would like to keep only Tp^ n _ 1 (D^ n _ 1 ), since the sum of these terms produce the 
drift we are looking for, but unfortunately at each step ( Z n ^) has the same order 

as Tp 1 p n _ 1 (D^ n _ 1 ). To get rid of this extra term p^, ri _ 1 (Z n ^), we use the fact that it is 
an increment of a martingale and thus averages to 0 under summation. More precisely, 
denoting z n := (Z n> TAT n ) and using Doob’s inequality we obtain, 

> n i/ 2 - 4 C e 


( sup 

J2 Zn 


1 ^ n ^ m 


z - 

1 ^ ri ^ rij 


(5.17) 


and we have the following decomposition: 
2\ 


E 


J2 z - 

1 ^ n ^ rif 


^E[ ]T E[|z n+ 1 | 2 |J- Tn 

1 ^ n ^ ny — 1 

tfC J2 E [\\ Z n,T^\W d ] < Cn f TN~\ 

1 ^ n ^ 72 y — 1 


(5.18) 




























28 


ERIC LUgON AND CHRISTOPHE POQUET 


where we have used dH3D- Since nj is of order N 1//2 , the probability in (15.171) tends to 
0 when IV —» oo and recalling (15.161) . we deduce that there exists an event 0^ satisfying 
P(^s>0 —>1 such that on n? 

n t 

= i/jq + tY^ Wn-iiD^n-i) + 0(A" 1/4+2C ) • (5.19) 


71—1 


The quantity ^ n -i( D ^n-i) = N 1/2 Pbn-i (-<% (6v • (J * tfVn-ikV’n-i)) depends linearly 
in and since the model is invariant by rotation, the projection does not depend on 
ip n _ i. So we can write it as N~ 1 ^ 2 b(^N), where the linear form b is given by 

HO : = p(-<% (€-( J *«)?)) = p ^ Y j • ( 5 - 2 °) 


We can rewrite (15.191) as 


'Z't* ~ Ho + 


T 


N 1 / 2 


N 1 / 2 


(t-$) 


6(6v) + 0(iV 1/4+2C ) • 


(5.21) 


Since 


t-t”- T 


u 0 n 1/2 




s 


AT 1 / 2 


and fo(£jv) = O(N^), we deduce 


= V>o + &&v)(* - O + 0(iV- 1/4+2C ), (5.22) 

which implies the result, since tjf = 0(N~ 1 ^ 2 logN). Proposition 15.21 is proved. □ 

We can now prove the following result, which together with Proposition 12.91 implies 
directly Theorem 12.31 

Proposition 5.3. There exists N sufficiently large such that, on the event ft!} , 


sup 

,tf] 


Tn,N 1 / 2 t ~ Qi> 0 +b(£ N )t 


—l,d 


= 0(lV" 1/4+2C ) 


(5.23) 


where the error 0{N 1 / 4+2 ?) is uniform on . 

Proof of Proposition 15.31 We place ourselves on the event introduced in the proof of 
Proposition 15.21 For each t such that A" 1//2 t € [T n ,T n+ i] we can decompose p N N 1 / 2 t as 

Tn,n^/H = Qti’n + u n+i,NVH-T„ ■ (5-24) 

But Proposition 15.11 implies that v n+l N i/i t _ Tn = 0(A r_1 / 2+2 ^) and for such time t we 
have 

(hl>n = = 9^, +Kt N )t + 0(iV" 1/4+2C ), (5.25) 

where we have used Proposition 15.21 □ 


6. Approaching the manifold 

The purpose of this section is to prove Proposition 12.91 We follow here the same ideas as 
in [7|, Section 5. From now on, we fix eq > 0 and po € Hf l such that dist^-i (po, M) ^ £q. 

a 

The parameter eo will be chosen sufficiently small in the following. We proceed in three 
steps: 

(1) We rely on the convergence in finite time of the empirical measure HN,t to the 
solution p t of (11.71) starting from po in order to show that fiN,t approaches M (up 
to a distance of order e 0 )- This step requires a time interval of order logeo- 
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(2) We use the linear stability of M under (11.71) and control the noise terms of the 
dynamics to show that the empirical measure approaches M up to a distance of 
order A _1 / 2 +2C. This step requires a time interval of order log IV. 

(3) We show that the empirical measure stays at distance IV _1 / 2+2 ^ from M up to the 
time ty . 

First step. As explained in Section 12.31 the stability of M implies that if Co is small 
enough the deterministic solution pt of the limit PDE (11.71) with initial condition po con¬ 
verges to a qe 0 € M. In particular, after a time si, pt satisfies ||p Sl — ( 7 e 0 1|—i,rf ^ £o- Due 
to the linear stability of M, this time si is of order — ^ logeo- 

In order to show that the empirical measure is close to the deterministic trajectory pt 
when N is large, we use a mild formulation similar to the one obtained in Section [2 but 
this time relying on the (2 d+ l)-dimensional Laplacian operator A^. More precisely using 
similar argument as in Section [HJ one can obtain the following equality in H^ 1 : 


VN,t~Pt = e2 Ad OuAT, 0 -po) - / e* L 

Jo 


do I Pn,s ® w + PN,t ^2 


* mL 


k=—d 


do lp s ®uj + p s ^2 xkj * P. 


k=—d 


ds + Zt, (6.1) 


where zt satisfies, for all test function / = (/ ,. • • > f d ) 


d \i pi 

= E w E [ d ° 

i=—d j =1 


/ 0 




K-(«))d B)(s). 


( 6 . 2 ) 


Since A^ is simply the classical one-dimensional Laplacian operator A on each coordinate, 
it is sectorial (in fact self-adjoint) with negative spectrum. Using the classical bound 


|e tA /||_i ^ -^=||/||_2 for the one-dimensional Laplacian operator, we directly obtain 


C 


< ^11/11-2,d, 


(6.3) 


and with similar estimates as the one used in Section [4] one can show that the event B^ 
defined as 


(6.4) 


Bi ■= { sup \\z t \\- 1)d < \ 


0 ^ t < si 


satisfies P(i?( v ) -A 1 as N -A oo. Let us write the shortcut 

U NyS . t := e q 6 jlNs (g) w ^ A k N J * n k Ns -dg lp s ®u + p s E ^ kj * Ps 


k=—d 


k=—d 


for the term within the integral in (16.11) . Note that the mapping (p, u) i-a dg{pJ * v) 
satisfies (see [7], Lemma A.3 for a proof) 


\\d e {p.J * u)||_ 2 ^ C'||/i||_i||u||_i 


(6.5) 
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Using (16.31) and ([6.51) . we obtain 


II d ~ 


e* Ad d e ( mat,. s E X x J *t4 , s ) ~ d e \Ps Y Xkj * p ‘ 


€ 


C 




c 


k=—d 

d 


k=—d 

d 


do ( PN,s Y A ^ J */4* p. E xkJ - 


k=—d 


k=—d 


d d 


+ 


c 


i=—d k=—d 

d d 


< = Y Xk de [Vn,s J * t4f,a) - d e (pU *P. 


-2 


— l,d 


—2,d 


\/t — S : 


E E * Mat, 


i=—d k=—d 






c 


-2 

c 


( 6 . 6 ) 

(6.7) 


- ui^a||-l,d + ||WV,s||-m)I|WV,s -Ps||-l,d + 1/2+C ||^JV,s||-l,(i 

( 6 . 8 ) 

(6.9) 


a 


MAI,. 


-PS ll-M + iV- 1/2+C ) , 


y/t — S 

where we have used in particular (12.41) . since both and //jv,s are probabilities. Let us 
place ourselves on the event 


: = 


| IIA*IV,0 — Po||-l,d < y } n B l 


( 6 . 10 ) 


which satisfies obviously P(iL^) -A 1 as N —>• oo. Then, for all t P si, (16.31) and (|6.6I) 
imply that (16.11) can be rewritten on the event as 

IImat,*-P t||-i,d ^ y + C \[^ N<: + c J 0 V T~ ~ ^ N ’ S ~^II-M ds > (6-H) 

so applying the Gronwall-Henry inequality ([24j. Lemma 7.1.1 and Exercise 1), one obtains 
that for some a > 0 (independent from N and £q), on the event B 2 and for all t ^ si 


11 Mm, t - Ptll-M ^ 2 (y + C \f^ NC ) ^ • (6-12) 

We deduce that for N large enough, the projection 

V>o : = Prej M (Miv,5i) 

is well defined and IImaAsi — 7hil|-i,d ^ £0 on B 2 , which means that |^q — #o| ^ C £ o and 
II Mai, si Q60 II ^ 2 e 0 - 

Second step. Now that we know that dist(/j,N S1 , M) ^ 2eq with increasing probability 
as N —» 00 , we can use a similar scheme as the one defined in Section T2.61 to show that the 
empirical measure approaches M up to a distance fV _1 / 2+2 ^ with high probability. Since 
this part is very similar to the work done in Section [5l we do not specify all the details. 

We consider the evolution of the dynamics on time intervals [T n , T n +i] with T n = si+nT 
where T is such that e~ 1IjT ^ 4 clc P • We consider also a sequence of real numbers h n 

satisfying h\ = 2 eq and h n + 1 = and take this time the number of step nj of our scheme 
as 

n/ := inf j n: h n ^ j\i^V 2 + 2 C j . 


(6.13) 
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It is clear that rif is of order O(loglV). To ensure the existence of the projections of the 
process on M at each step, we introduce, as in Section [2761 the stopping time 

(n T ,r) = inf{(n,f) <E {1,... ,n f } x [0,T] : \\Hf n _ 1+t ~ 9a n _i ||-i,d ^ c} , (6.14) 

where a n = proj M (^ ) when it exists. This allows us to define the random phases ipn-i 
defined as 

:= proj M (Ai (nASr _ 1)f+tA ^), (6.15) 

and the processes v n ^ dehned for n = 1 ,,rif as 

v n ,t ■= M( nA ^ T _i)f + t A ^n - • (6.16) 

This last process satisfies the mild equation 

^ n ~ 1 (D^ n _ i + R^ n _ i (i^ ritS )) ds + Z n ^n , (6.17) 


where Z Ujt is defined as 

_ d i Ni r t 

Z n,U) = E wT. a » 

i=—d j=l J0 

Section [I] shows that the event 

A N = < sup ~ sup^ \\Z njt \\-i,d < f l/2 N~ l/2+c - \ , (6.19) 

satisfies P(A iV ) -» 1 as IV —>• oo. 

In the first step of this proof we have shown, since tj )o = i/Iq, that, on the event 
(recall (16.161) 1. we have ||vi,o||-i,d = ||/Cv,si — <fyi||-i,d ^ h\. Our aim is to prove that on 
the event B% defined as 

B$ := A N HB? , (6.20) 

we have ||un,o||-i,rf ^ h n for all n = l,...,rif. This would imply, using the notations 
s 2 = T nf and i/’o = proj M (///v, S2 ), that \\J1n,s 2 ~ Q^\\-i,d < IV _1/2+2C . We place ourselves 
on the event B^ . From the mild formulation (16.171) . if n < rif and H^oII-m ^ we get 

||^||-M < C L e-^ n h n 

+ 2C L T\\Dr || — i,d + C L (t + 2T 1 / 2 ') sup ||i?r ( v n ,s)\\-i,d 

vn v ' 0 < s ^ t 

+ r 1/2 IV” 1/2+c . (6.21) 

Consider the time t* defined as 


(t—s)L*~ 

e / 


(^(T n _! + a)) d^(T n _! + a). ( 6 . 18 ) 


■'n.t 


= e 


[t/\T n )L 


V’ra — 1 // n — 
u n,0 


rt/\T rl 

Jo 
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For all t ^ t* we have 

SUP 1 (* 7 n, S )||-l,d 

0 sC s < i 

^ Cl SUp \\Vn,s\\-l d + N ~ 1/2 max 1^1 SU P ll^n.sll-M 
V 0 < s A t ’ k=-d,...,d o^s^t 

< C{C 2 L hl + C L N~ l ' 2+ th n ). (6.23) 

The last quantity is smaller than C(N,So)h n , where C(N,£q) -A 0 as N —>• oo and £o —> 0. 
On the other hand we have shown in (15.31) that \D-~ i ||_i ! d ^ CN 1 / 2 +?. Since n <rif 

we have hjsr > 2 jJr-IV _1 / 2+2 £, which means that IV -1 / 2 " 1 "^ is negligible with respect to h n for 
N large enough. So for N large enough, t* ^ T and we have (recall that e~ XT ^ 4 clc P ) 

1 3 

II^tII-M ^ ACp h n + ^ s Cp h n ’ (6.24) 

when £0 is small enough. It remains to show that ||Pn+i,o||-i,d ^ to conclude the 
recursion. We do not prove it in details, since it can be done by proceeding exactly as in 
the proof of Proposition 15.11 decomposing \\i/ ^||_i,d and showing that it can be written 
as 

ll^ra+poll — l,d ^ M t (6.25) 

which implies that |]Pn+i,o||-M ^ f^n + 0 (h 2 ) ^ on the event when £0 is small 
enough and concludes the recursion. Note that the estimate for ip n — ip n -i obtained in 
452D leads to 

n f n f 

IV’o-V'ol < J2 \'ipn 'ipn — 1 | ^ C ^ ^ hn ^ 2Ch\ ^ A.C so , (6.26) 

n= 1 n =1 

on the event which gives | -0o ~ $o| ^ C'e 0 for some C'. 

Third step. In the previous step, we have constructed a time s 2 such that s 2 ^ — 
j log £0 + Ci log N for some constant Ci and such that ||jujv,s 2 — ||-i,d ^ N _1 / 2 + 2 C with 

high probability. We can now consider a time s 3 = clogIV for c = Ci + 1, which does 
not depend in £q. For N large enough, we obviously have S 3 > S 2 - In order to prove that 
ll/Mr,s 3 - ll-M ^ IV _1/2+2C with high probability, where = proj M (/xjv,s 3 ), it suffices 
to decompose the dynamics on the interval [S2, S3] according to an iterative scheme with 
time step T satisfying e~ lLT ^ 4 ClC p as does T and apply exactly the same procedure as 
in Proposition 15.11 

This last step induces a phase shift \ipQ — ?/>q| ^ CN~ 1 / 2+2 C l G g N ^ Ce 0 , for N large 
enough. This concludes the proof, with = IV _ 1 / 2 S 3 . 

7. Estimates on the drift b 

7.1. The case of a symmetric disorder. We prove here Proposition 12.41 and drop 
for simplicity the dependency in ip and 5. We consider £ = (£ _<i ,..., £ d ) such that 
£ _ * = £* for all i = 1 ,d and aim at proving that &(£) = 0 , where the drift &(£) = 

p (~ d e ({EL-d^ fc ( J * 9 fc )} 1 )) is g iven b y »• 

The space of regular (C 2 , say) test functions / = (f~ d } f~( d ~ 1 \...,f d ~ 1 ,f d ) can be 
naturally decomposed into the direct sum of the space O (resp. 8) of odd (resp. even) test 
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function in both variables (9, i), that is / E O (resp. / € £) if and only if f~ l (—9) = —/*(#) 
(resp. f~ l (—9) = /*(#)) for all 9 E T and i = 0,..., d. One easily sees from the definition of 
</(■) in (11.41) and the definition of g in (| 1.8 1) that q E £ and ((J * q ~ d ),..., (J * q d )) E O. 
Let us denote Q{9) := Ylt=-d£ k (J * Q k )(0)- Using that = £*, one obtains that 
Q(6) = £°(J * q°)(9) + J2k=i^ k ((J * V k )(+ (J * Q~ k )(9)), so that we deduce that Q is 
an odd function of 9 and that 9 i-A Q(9)q(9) E cO. Consequently do {Q{9)q{9)) € £. 
Hence, in order to prove Proposition 12.41 it suffices to prove that 

Wi E E, p (h) = 0. (7.1) 

This is indeed the case since one easily sees from the definition (12.131) of the operator 
L = L^, t s that L(S) C £ and L(0) C O and since p is the projection on the eigenfunction 
doq E O. Proposition 12.41 is proved. 

7.2. Small 5 asymptotics of the drift. Our aim here is to prove Proposition 12.51 that 
gives the first order expansion of the drift f>(£) defined in (|5.20p as 5 —» 0. Due to the 
rotational invariance of the system, we can work with the stationary solution go ,5 that we 
denote q$ throughout this section. We denote p^ as p^=o ,<5 (recall ([2.151) 1 and D$(£) as 
Dn, o, 5 , (recall (12.22 [) ). With these notations the drift b is given by 

b(0 = P s(D s (0). (7.2) 

When 5 = 0, it is straightforward to see that qs = (qj d , ■ ■ ■, q d ) is equal to (go, • • •, go), 
where go is the stationary solution of the nonlinear Fokker-Planck equation without dis¬ 
order (12.121) . We refer to Section lB.ll below for precise definitions (see in particular (IB.2I) 
and (IB.3D where the normalisation factor Zq and the fixed-point parameter ro are defined). 
The following result (proved in Appendix[D|) provides the next order of the approximation 
of qs as 5 —>• 0. 


Lemma 7.1. For i = —d ,..., d we have 

ql(9) = q 0 (e)+Su i K(9)q 0 (9) + O(S 2 ), 

where 

f„ 27r e~ 2Kro cos “ d u f n 2n e 2Kro cos U udu 

k(9) = 29 + 47t—---24°- 


(7.3) 


Zn 


r 2?r e 2Kro cosrt r 2?r e ~2Krg cos v ^ 

- 47t—--, (7.4) 


Z 2 

"0 


and where the error 0(5 2 ) is uniform in 9 G T. 


The projection p^ also converges in some sense to the projection po on the tangent space 
of the stable circle of stationary profiles of (12.121) at go- Moreover, the system given by 
(12.121) admits a nice Hilbertian structure, which allows to know po explicitly. This allows 
us to obtain the following first order expansion of p^, whose proof is given in Appendix [Dj 

Lemma 7.2. For all coordinate by coordinate primitive (U~ d ,... ,U d ) of u smooth, we 
have 

-72 d 


P«(«) = 


-'0 


Zl - 4vr 2 

u k=-d 


E*‘/ 

j 4T 


1 - 


2t r 
Z^qo 


U k + 0{5\\u\U 4 ). 


(7.5) 


We have now the tools required to obtain the first order expansion of the drift &(£). 
The result we want to prove is 
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Proposition 7.3. For all £ such that Ylt=-d^ k = 0 we have 


b(0 = 8 J2 Z k u k + 0(6 2 ). 


(7.6) 


k=—d 


Proof of Proposition \7,3\ First remark that when <5 = 0, we obtain, using Lemma f7.ll that 
for all i = —d,..., d: 

d 


Dh(0 = d e 


<7o 




* <7o 


k=—d 


= o, 


(7.7) 


since Y2t=-d^ k = 0- We deduce, using again Lemma [7.11 the following expansion for 

Dm 


D}(0 = 8u l d 0 


nq 0 


E 


* <70 


k=—d 


+ 8d e 


<7o 


E m k J*(*q o) 


k=—d 


+ 0 ( 8 2 


= 8d e 


<7o 


£ m k J*(*qo) 


k=—d 


+ 0(8 2 ), (7.8) 


where we have used again the fact that Ylk=-d^ k = 6- Applying Lemma fL2l we deduce 

d d 


HO = <5 


z 2 


L Z 2 - 4vr 2 


1 - 


2vr 

Zfao 


q 0 J * (k<?o) 


E E *m k + 0(8 2 ), (7.9) 


J i=—dk=—d 


and recalling that Yj=-d = 1 and denoting 


c b : = 


7 2 


— 47r 2 


/ 
J T 


1 - 


2vr 

Zoqo 


q 0 J * (K,q 0 ), 


we simply obtain 


6(0 = Sc b J2 m + O(8 2 ). 


(7.10) 


(7.11) 


k=—d 


It remains to show that c b = 1. Now using the fact that J(9 — 9') = — IF sin# cos 0' + 


K cos 9 sin#', f 2w sin(0)<7o(O = 0 and cos(9)qo(9) = vq, we obtain 


2tt 


and 


/•27T /*27T 

/ q 0 (9)J* (Kq 0 )(9)d9 = Kr 0 / sin(9')K(9')q 0 (9') d9', 

Jo Jo 

r2i r 

/ J * (nqo)(9) d9 = 0 . 

Jo 


So the constant c b can be simplified as follows 


c b = 


Kr 0 Z , 
Zq — An 2 


2 r 2tt 


which leads to 


c b = 


2Kr 0 Z 2 
Zq — 47t 2 


L 


27T ^2Kvq cos 0 

sin C 


/ sin(9)K(9)qo(9) d9 , 

Jo 

fg™ e~ 2Kr ° cosu du 


Zn 


9 + 2tt- 




d 9 


(7.12) 

(7.13) 

(7.14) 

(7.15) 
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Integrating by parts and using the fact that dg[e 2Krocos ^ 9 ' > } = —2Krosin9e 2Kr ° cos9 , we 
obtain 

/•2tt 2R>0 jr 

/ 9 sin 9e 2Kr ° cos 9 d9 = ---+ —^ , (7.16) 

Jo Kr 0 2 Kr 0 K J 

and 

f 2 ” sin 9e 2Krocos 9 e- 2Krocosu du = - -L- , (7.17) 

Jo Je 2Kr o K r 0 

which implies that (^27777 — oHrlfz 7 ) = * • Proposition 17.31 is proved. □ 

Using Proposition El we can now compute the first order of the variance v 2 of the 
limiting normal distribution of 6 (£/v) when the disorder is i.i.d: 


Proof of Proposition \2.5\ From the Central Limit Theorem, we know that £at converges 
as IV ^ 00 to a Gaussian distribution with mean 0 and covariance matrix £ satisfying 


r £ fe , fc = A fc (l — A fc ) ke{-d,...,d}, 

\ E fc ,{ = —X k X l k,le{-d,...,d},k^l. 

Applying Proposition 17.31 we obtain 


(7.18) 


v 2 = 5 2 | A fc (l-A fc )(, 

. kE{—d,...,d} 


c/) 2 - 


Y X k X l uj k u l ) + 0(6 3 ), (7.19) 

k,l£{-d,...,d},k^l 


and since A k = X k and u k = —u k the terms with l 7 ! —k cancel in the second sum, 
which gives the result. 

□ 


Appendix A. Construction of rigged-spaces 

We specify here the construction of the Hilbert distributions spaces we work with in 
this paper. It is based on the notion of rigged Hilbert spaces (see Ha,?, si). 


A.l. Functional spaces on T. Consider Lq := {-a € L 2 , J T u{9) d0 = 0}, the space 
of square integrable functions with zero mean value, endowed with the norm ||rt || 2 := 

(f T u(9) 2 d9) 5 . We call a weight any strictly positive function 9 i-a w(9) on T. For any 
weight w on T, define H ^ as the closure of {u £ C 1 (T), f T u(9) d 9 = 0} w.r.t. the norm 

\H\i, w ■= ( d e u(9)) 2 w{9) d6^ . 


There is a continuous and dense injection of into Lq and the corresponding dual space 
can be identified as H~,, that is the closure of {u € C 1 (T), Jju(9)d9 = 0} under the 
norm 


\u\ 


— 1,1 /w 


m 2 

w(9) 


d 9 


where IA is the primitive of u such that f T ^ = 0. 
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A.2. Functional spaces on TxM. The correct set-up of the paper is to consider test 
functions of both oscillators and frequencies, that is (0,uj) e -)- u(9,ijo), where 6 £ T and 
uj € R. Since the disorder is assumed to take a finite number of values ... , u; rf }, 

it is equivalent to consider vector-valued test functions 6 i-a (u~ d (9),... ,u d (6)) and it is 
straightforward to define the counterparts of the norms defined in the last paragraph for 
these vector-valued functions: Consider Lq^ := (Lq) cZ endowed with the product norm 




In the same way, consider the space d , closure of {(u d ,..., u d ) £ C 1 (T), f T u k {6) d0 = 0} 
under the norm 

i 

/ d 


I l.w.d 


E A ‘ 


\k=—d 


U 


as well as the space H 1 J w d endowed with the 


norm 


u 


— 1,1 /w,d *" 


E A ‘ 

\k=—d 


U 


2 

1 ,w 




(A.l) 


(A.2) 


Note that if w\ and W 2 are bounded weights, the norms l]-^ and ||-1| j W2 (resp. ||-||_ 
and ||• || i/ W2 ) are equivalent. The same holds for the (2 d + l)-dimensional norms. 


A. 3. Fractional spaces. Define also the fractional norm || • || a ^ (where a ^ 0): consider 
Arf the Laplacian operator on each coordinate, || • ||o the L 2 -norm on T and ||u||orf = 
E fc A fc ||u fe || 2 and define 

IHIm = = E VII(1-A)“/V||J. (A.3) 

k=—d 

We denote as H d the closure of regular functions with zero mean-value on T under the 
previous norm and H d a the corresponding dual space. 


Appendix B. Spectral estimates and regularity results on semigroups 

The purpose of this paragraph is to establish spectral estimates on and its adjoint 
as well as regularity estimates on their semigroups e and e 

B.l. The case 5 = 0. The analysis of the dynamics of (11.51) and (11.71) is based on pertur¬ 
bations argument on the mean-field plane rotators system (12.111) and (12.121) . The proof 
relies in particular strongly on the fact that (12.111) is reversible, with an explicit free en¬ 
ergy ISH5I- However, one should note that the limit as 5 —» 0 of (ESI) or E3> is slightly 
different to the mean-field model (12.111) - (12.121) . In particular, (11.71) becomes as 5 —>■ 0 

dtql(0) = ^d 2 e q\{0)-d e 5Z P kj *9tW^ , i = -d,...,d, (B.l) 

which corresponds to the situation where the disorder is no longer present but where 
the rotators have been (artificially) separated in different subpopulations. Following the 
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terminology of m where (IB.ID has been already encountered, we call this system the non- 
disordered system. It is shown in [21] , Section 2.1, that the non-disordered system (jB.lD 
presents most of the properties of the mean field plane rotators model (|2.12l) . In particular, 
for all K > 1, one can show that (1B.1|) admits a unique circle Mo t nd of synchronized profiles, 
that is stable as t —>• oo. Mq is given by the translations of the profile qo,nd = (qo, • • •, qo), 
where qo is the profile generating the stable circle Mo of non trivial solutions of ()2.12|i . 
namely 


g2Kro cos 6 

q ^ 9) := Zo(2Kro) 


(B.2) 


where Zo(x) = 2 ttIo(x), Iq(x ) = ^ Jq™ e xcos ^ d 9 is the standard modified Bessel function 
of order 0 and tq is the unique positive solution of the fixed-point problem 


r 0 = ^o{2Kro ), with To(x) 


cos(6)e xcos 9 d0 
Z 0 (x) 


(B.3) 


The derivation of these stationary solutions is highly similar to the procedure described 
in Section [131 and we refer to the aforementioned references for more details. Note that 
one can draw a simple correspondance between the present definitions and the definitions 
of Section [131 in the case of 5 = 0: namely, one readily sees that, for any i = —d ,..., d, 
Sq(9,x) = e xcos< '°' 1 Zq(x) (recall (jl.9[> ) and Zq(x) = Zo(x) 2 , so that the definition of 
when 5 = 0 (recall (II. IIP ) coincides with Tq given in (IB.3D . 


B.2. Spectral estimates when 5 = 0. Define the linearized operator around any sta¬ 
tionary solution q 0 ,nd £ M 0 , n d- 


(•AuY = Ijdjju 1 - dg f (J * q 0 )u l + q 0 ^ J * u k j , i = -d ,..., d , (B.4) 

with domain V(A) = {(u~ d ,... , u d ) G C 2 (T) 2d+1 , f T u k (0) &6 = 0, k = — d,... , d}. We 
recall the following result (see [21], Proposition 2.1): 


Proposition B.l. A is essentially self-adjoint with compact resolvent in P/) , d . Its 
spectrum lies in (—oo,0], 0 is a simple eigenvalue, with eigenspace spanned by dgqo, n d- 
The spectral gap between 0 and the rest of the spectrum is denoted as 7 a ■ 


One can deduce from Proposition IB. 1 1 similar spectral properties of its dual A* in L 


2 . 
0 ,d‘ 


(■ A*vy := ^dgV 1 + (J * qo)dgv l - f {{J * qo)dgv l ) dd-^2 X k J * (q 0 d e v k ), i = -d ,..., d , 

k =1 

(B.5) 


with domain V(A*) = T>(A). 


Proposition B.2. A* is essentially self-adjoint with compact resolvent in H\/ qi> g- Its 
spectrum lies in (—00,0], and 0 is a simple eigenvalue and its spectral gap 7 a* is equal to 
7 A- 

Proof of Proposition \B.2[ Let us introduce the operator U defined from H^. d to Hf} qo g 
as 


Uf(9) := -d e (q 0 (6)d e m). 
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U is an isometry between d and U realizes a bijection from {u € C°°(T) rf , 

f T u k (9) dd = 0, k = 1,... , d} into itself and for every /, g G H qo d , 

/rr , TT , v- f h(o)def k (0)) (, qo (e)d e g k (e )) 

= ?/ T -iw- 


dfl 


= E / ?o(0)^/ fc (0)^ fc (0)d0 = (/,</> 1>m>d • (B.6) 

fc 

Moreover, the following identity holds: 

A* = U~ l AU , (B.7) 

so the operators 7l on d and A* on H\/ qn d have the same structural and spectral 

properties. □ 

B.3. Spectral estimates of $ and its adjoint. We are in position to deduce spectral 
estimates on the disordered operators L$ and its adjoint L* s in Lq d (we drop the index if 
in this section for simplicity). 

Proposition B.3. The adjoint L* d of L$ in Lq d is given by for all i = 1, ... , d 


1 “ 

(L* sV y = d^ + S^de^ + ide^) Y A fe J*g£ 


k=—d 


(de^) Y xkj *q k ) de 

s. k=—d ) 

d 

-Y X kj *(q k 5 d ev k ), (B.8) 


k=—d 


with domain D(L* d ) = D(A). 

Proof of Provosition \B.SX For all regular u and v, 

{L* s v,u) 2 ,d = (v,L 5 u) 2 ,d 


Y A* ( v\ \d 2 e u* - SJdou* -dei^Y xkj * qs +qsJ2 Xkj * 


u 


i=—d 
d 


k=—d 


2 ,d 


= Y xi ( \ d e yi + Sujid 9 yi + d o vi Y xkj *1s- l I 9<)V l Y J 

i=-d \ 2 k=-d M 

d d 

+ Y Y A 4 A k (qid e v\J*u k ) 24 

i=—d k=—d 

= /^P + S^deP + ideP) Y xk J*Q k ~ [ ((dev*) Y 

\ k=—d V k=-i 

~Iy a k J*(q k d e v k ),u^) , 


k=—d 
d 

*qsjdd,u l 

k=—d / / 2 ,d 


A k J*q k 5 dd, u* 


2 ,d 


\k=-d / 2 ,d 

which precisely gives (13.51) . 

The main result of this section is the following 


□ 
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Proposition B.4. There exists 62 = ^(K) > 0 such that for all 6 ^ 82 , everything that 
follows is true: the operator L * 5 (resp. L g ) is sectorial in H qo d (resp. d ), its spectrum 

lies in a sector of the type {A E C : | arg(A)| > 7 r /2 + a} for some a > 0 and 0 is an 
isolated eigenvalue for L* d (resp. Ls), at a distance from the rest of the spectrum denoted 
by 7 l* (resp. 7 l s )■ Moreover, both Ls and L * 5 generate a Co-semigroup t i-a e tLs (resp. 
t i-a e tL s)in Lq d and e tL s = (e tLs )*. 

Proof of Proposition \B-4\ The result concerning the operator Ls has been proved in [2Tj, 
Th. 2.5. For the sake of completeness, we recall here the main arguments concerning 
in H qo d but we refer to [21], Section 6.2 for precise details. Note that we need a precise 
control of the spectrum of L d around the origin. In particular, one has to ensure that 
the spectrum of L* s remains in the negative part of the complex plane. We write L* d as a 
perturbation for small disorder of the non-disordered case: 

L* S = A* + B S , (B.9) 

where A* is given in (IB. 51) and Bs is a small perturbation as 8 —>• 0. More precisely, 
following the exact same strategy as in [21], Proposition 6.5, p. 356, one obtains that the 
operator Bs is A*-bounded: there exist constants as and bs (only depending on 8 and K) 
such that for all u in the domain of (the closure of) A* 

W B M\i , qo4 < a s\\ u \\l, qo ,d + b s\\ A * u \\l, qo ,d’ ( B - 10 ) 

with as = 0(8) and bs = 0(5), as 8 —>■ 0. Note that the only things that differs between 
this result and [21] . Proposition 6.5 is that we work here with an Pf 1 -norm whereas the 
result in [21] concerns an iJ -1 -norm. 

Fix some e > 0 (that will be specified later) and define L* s := — e and A e := A — e, 

so that Lg e = A* + Bs■ Fix a E (0, ^) and introduce the following subset of the complex 
plane 

:= |A E C, |arg(A)| < ^ + a j U {0}. 

The operator A e (as A itself) is self-adjoint in H_ 11 / qo and hence, sectorial. In particular, 
there exists M > 0 such that ||i?(A, A e )\\ H ~i ^ pq, for all A E S Q . Note that the 

constant M is indeed independent of e > 0 and that the previous inequality is also true 
for A in place of A e (see [21], (6.12)). Using (IB. 71) . one obtains that ||i?(A, A *)|| h i ^ ttt. 

90 , d 

For A E u E H l qo d , 

\\ B sR(X, A*)u\\ l qo d ^ as ||7?(A, A*)u\\ lqod + bs ||A*i2(A, A*)ti|| liJOj(i , 

^ ~\>f + ( M + Nil, 90 ,d ’ 

Choose 8 sufficiently small so that b$(l + M) ^ | and ^ j. Then for |A| > e ^ 4 Mas, 
we have \\BsR(X, A*)u \\ 1 ^ f x so that the operator 1 — B g R(X, A*) is invertible from 
Hq od into itself, with norm smaller than 2. A simple computation shows that in this case 

(A _ (A* + Bs ))- 1 = R( A, A*)(l - B 5 R( A, A*)) -1 , 
which gives that, for A E £ a , |A| > e, ||i?(A, L* s )\\ H i ^ 44f. Consequently, the spectrum 
of Lg is contained in 

©a, £ := |aeC, ^ ^ arg(A) ^ ~ « j 


U{AeC, |A| ^e}. 










40 


ERIC LUgON AND CHRISTOPHE POQUET 


In particular, 0 € p(L| 2e ) and for all A € C with 3J(A) > 0 (hence |A| < |A + 2e|), 

R(KLUJ 


H 1 , 
QO 




5,2e- 
M 

]A+2e| 




. The fact that this estimate can be extended to some E c 


for some a! is a consequence of a Taylor’s expansion argument (see m , Proposition 6.2), 
so that L* S2e (and L*f) is indeed sectorial. 

At this point, we cannot rule out the possibility that some elements of the spectrum of 
L* s may lie in 0 £iQ ,n{A G C, 5R(A) > 0}. The last point of the proof is to show that one can 
choose e and a smaller d such that this situation does not hold: choose e = ^ > 0, where 
7,4 is the spectral gap of A. In particular, the circle centered in 0 with radius e separates 
the eigenvalue 0 (of multiplicity 1) from the rest of the spectrum of A*. An application 
of [26], Theorem IV-3.18, p. 214, shows that one can choose 5 sufficiently small so that 
the spectrum of the perturbed operator is likewise separated by this circle: for such 6, 
there is a unique eigenvalue (with multiplicity 1) within the boundary of this circle). But 
we know already that 0 is an eigenvalue for the perturbed operator L* s . By uniqueness, we 
conclude that there is no eigenvalue in the positive part of the complex plane. We leave 
the details of this argument to [21], Section 6.2.5. 

Using [36], Corollary 10.6, p. 41, L* s is the generator of the adjoint of t e+ e tLs in Lq^, 
which is a Co-semigroup. This concludes the proof of Proposition IB. 41 □ 


B.4. Equivalence of norms. For any 0 ^ /3 ^ 1, consider the interpolation norm || • \\yp 
associated to the sectorial operator 1 — L* s defined as 

My? = \\(l-L* 5 fu\\ liqo>d . (B.ll) 

Recall also the definition of the fractional norm in (IA.3I) . 


Lemma B.5. Under the assumptions of Proposition B. 4 , for any 0 ^ (3 ^ 1, there exists 
ci, Ci >0 such that for all u, 

Cl|M|i + 2/3,(2 ^ IMIv/ 3 ^ Cl IMI1+2/3,d ■ (B-12) 

Proof of Lemma \B.5[ We can decompose L* & as follows: 


L* 5 = -A d + R , 


(B.13) 


where, for all i = 1,..., d 


k=—d 


k=—d 


(.Rvf = w + v J2 \ k J*q k Q - J2 A k j*{q%{d e v k )) 

dev^e) ^2 x k j*q k (e)] de. (B. 14 ) 


k=—d 


Since R only contains first order derivatives and J and Qq are smooth, it is easy to see 
that for all u G H%, we have 

II-R^IIm ^ C'||tt|| 2 ,d- (B.15) 

One deduces immediately from this estimate that there exists a constant C > 0 such that 
for u £ Hj, 

|| [2(1 — Lg) — (1 — Ad)] u|| 1)d <<7|M| 2id . (B.16) 

Consequently, the operator [2(1 — L* s ) — (1 — A^)](l — A^) -1 / 2 is bounded in H\. Since 
1 — Lfi is sectorial in H J with the same domain as Ad, an application of [21], Theorem 
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1.4.8 shows that the norms ||(1 — d and ||(1 — d are equivalent. The norm 

equivalence (IB.121) follows directly from the definitions (IB.Ill) and (IA.3D . □ 

B.5. Regularity of semigroups. Recall here the definition of the projection P° s on the 
kernel Span (dgq^j) of L^ ,5 dehned in Section [231 We drop here the dependance on if for 
simplicity. The corresponding projection on the kernel of Lg is given by P?’*. This kernel 
is one-dimensional, spanned by some 9 1 —> vq{ 9 ) and there exists a linear form p, bounded 
on H d such that, for all u E H d , Pg’*u = p(u)v o- Note that it is easy to see that vq is a 
regular (C°°) function on T. 


Proposition B.6. Suppose the assumptions of Proposition B.f are true. For any 7 E 
[0, 1l*), & n y P € [0,1] and all t > 0, u € H d , 

e-H 


tLZ / -1 nll.tx 

e *(1 -P s ’ )u 


. 0 ,* > 


1+2/3, d tP 


7-)0,*\ II 

- P 5 M|l,d, 


and 


t T * 

e tL s u 


1+2/3, d 


<c[l + 


W 




(B.17) 

(B.18) 

(B.19) 


and for all (5 ^ 0, ft' ^ 0 such that /3 + ft' ^ 1 and all h € p[^ +2 9+ 2 h , 
e tL s -ljh d =+ ^11(1 - P,5 0 ’*)h||i + 2/3'+2/3,d ■ 

Proof of Provosition lK6 1 Following Proposition IB.41 LgPg’* = 0 and L* s ( 1 — Pg’*) is sec¬ 
torial in H d , with spectrum lying in {A € C : | arg(A)| > ir/2 + s'} — 7 l* for some s' > 0. 

Let us first prove (IB.171) and (IB. 181) . Using [21], Theorem 1.4.3, (recall that 7 < 7 l|)j 
we obtain that for all t > 0 : 

ll ( _ L *yyu;(i_po,*) n || M ^ C '^ e -^ l |«|| lid . (B. 20 ) 

Now as in the proof of Lemma IB. 51 we can apply [21], Theorem 1.4.8 to show that the 
norms induced by (—Lg)P and (1 — Lg)@ are equivalent on the range of (1 — Pg’*) and we 
obtain for all u € H d 

\\e tL ht\\ l+2M < C\\{l-L* 5 fe tL HP 0 5 '*u+{l-P° 5 '*)u)\\ 14 < C^(l + e -^)||«|| M . 

(B.21) 

We have used here in particular the fact that for all u E H d , 


e tL sPg’*u 


1+2/3, d 


< I p(u) 


t j * 

e * Vo 


1+2/3, d 


= p[u 


since H't'oIIi+ 2/9 < + 00 . Concerning (IB.19I) . remark that 


e tL s - 1 = ( e tL, ~ s (1 - Pf*) - 1 ) (1 - P 


JL 


0,*N 


1+2/3, d ^ C IMIl,d > 


?■*), 


(B.22) 


so applying Theorem 1.8.4 of [21] we have 


e tL * s - 1) h 


1+2/3', d 


< C 


1 ,d 


(1 - L' s f (e“<*(l - Pp) - l) (1 - Pp)h 

<- C.y ||(1 - - Pp)h || M < Cpf >||(1 - Pp)h\\ 1+2lI , +2M . (B.23) 

This concludes the proof of Proposition IB.61 □ 

One can deduce from Proposition IB. 6 I a similar regularity result concerning the semi¬ 
group t e+ e tLs : 
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Proposition B.7. For all K ^ l, all 0 ^ 5 < 5(K), the semigroup 1 e+ e tLs is continuous 
from Hf 2 to : for all h £ Hf 2 , t > 0, 

IK L ^||_ m <c(l + -j=y\h\\_ 2jd , (B.24) 

and for all £ £ ( 0,1/2), t > 0, u ^ 0, 


e (t+u)L s h _ e tL s h 


-1 ,d 


st Cu £ 



tV 2 + £ 



(B.25) 


Proof of Proposition \BT7 Let j5 £ [0,1], t > 0, h £ H d 1 and v a regular test function. 
Consider (hi)i ^ i a sequence of elements of d converging to h in . For all l ^ 1, 


\(e tL 'hi, v) d \ = {e tLs h u v) 2i 
I-(1+2/3),d 


hi , 


2,d 


€ 




1+2/3,d ^ C ^ hl Ul+2P), d ( 1 + Jh) Hi 


d ! 


where we used (IB. 18[) in the last inequality. Since hi converges to h in if -1 , one can make 
l -+ oo in the previous inequality and obtain |(e tL,5 /i, ^ C ||/i||_/ 1+2j g) ^ (l + jj?) d 

and since this is true for all regular v, one deduces that 


=“‘'‘11-1., «o(i + l) 


I-(1+2/3), d- 


(B.26) 


which gives (|B . 1 7D when /3 = i. In the same way, an immediate corollary of (IB.191) is that 
for all (5 ^ 0, (3' ^ 0 such that (5 + (3' ^ 1, for all t > 0 


\{e tLf - 1) h\\_ 


(l+2/3+2/3'),d ^ ^ II"'ll-(1+2,0 '),d 


We now turn to the proof of (IB.251) . Fix s £ (0,1/2) and apply (IB.26D for /3 
and (IB.271) for /3 = £ and fF = ^. 


(B.27) 
1/2 + e 


e (t+u)L s ^ _ e tLsfo 


-i,d^ C V 1 + t 1/2+£ 

^ Cu e ( 1 + 


o uL S 


-1) All- 


(2+2e),d : 


t l / 2+e 


1 - 2 ,d ■ 


This concludes the proof of Proposition IB. 71 


□ 


Appendix C. Projections 

The purpose of this section is to prove several regularity results concerning the pro¬ 
jection P° 5 u = p ii),s{ u )deq ^,5 (recall Section [2731 and ([2.151) 1 and the projection on the 
manifold M proj M (-) defined in Lemma 12.81 

Proof of Lemma 13. SI We first prove that ip e+ p^, is smooth. This follows from the fact 
that the whole operator is regular in ip £ T: we prove indeed that the mapping ip H > L^ 
is in fact real holomorphic, in the sense of Kato [26], p.375. Since the problem is invariant 
by rotation, it suffices to study the regularity of is a neighborhood of ip = 0. From 
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the definition of the stationary solution q in (11.81) , it is straightforward to see that one can 
expand in series of ip around '0 = 0: 

= q o (0) + JT d ^U =0 (^- 

k^l 

From this expansion, one deduces a similar expansion for around ip = 0: for all / 
regular 

WJ = L 0 f+J2 ^ k U k f , 

k > 1 

where each Uk is a differential operator of order 1, so that each Uk is relatively-bounded 
w.r.t Lq. In particular the hypotheses of [26], Theorem 2.6, p. 377 are satisfied. In 
particular, (L^)^ forms a real-holomorphic family. In particular, the mapping ip i—>- P® 
is also regular m, Theorem 1.7, p. 368), and so is the mapping ip i-a p^,. Then the 
mapping f(ip,h) = p^(h-q^) satisfies for each fixed ip 0 , f(ip 0 , Qp 0 ) = 0 and d^f{ip 0 , q^ 0 ) = 
—Vi’o^i’qtpo = — 1- So by the implicit function theorem, for all h in a certain neighborhood 
of gy, 0 , there exists a unique ip =: proj^/i) such that f(ip,h ) = 0 and h i->- pio] M (h) is 
smooth. □ 

The next result states that the first order of the projection proj M around q^ is given 
by the linear form p^, defined in (12.151) . 

Lemma C.l. For ip € T ; h € Hff 1 such that proj M ((jy, + h) is well-defined, we have 

proj m (QiI> + h) = ip + p y(h) + 0(\\h\\\ d ) . (C.l) 

Proof of Lemma \C. 1[ Consider the real u such that proj + h) = ip + u. Due to the 
smoothness of proj M , we have u = 0(||^i||_i,d)- The real number u satisfies 



(C.2) 

A first order expansion leads to 


p ^{h-ud^q^) = 0(u 2 ), 

(C.3) 

which gives the result, since p ip(d^q^) = 1. 

□ 


Appendix D. Expansions in <5 

The aim of this section is to obtain first order asymptotic of the drift in Theorem 12.31 
for small 5. We use the notations qg, p^ as in Section 17.21 putting the emphasis on the 
dependency of the different terms in 5. We denote also as rg > 0 the unique positive 
solution to the fixed point relation rg = ^g(2Krg) (recall ( 11 . 101 ) ). We begin with a result 
concerning rg as 5 —» 0: 

Lemma D.l. The mapping 5 >—)• rg is C°° and its derivative r'( 0) at 8 = 0 is zero, so 
that as 5 ^ 0: 

r s = r 0 + O(8 2 ), (D.l) 

where ro is the unique non-trivial solution of the fixed-point problem without disorder 
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Proof of Lemma \D.l[ Consider the C°° mapping g(r,5) = ty$(2Kr) — r. This mapping 
satisfies d r g{r\). 0) = 2Kd x ^o(2Kro) — 1. The fixed-point function r i-a To(2A>o) is strictly 
convex when AT > i m, Lemma 4), with derivative at the origin strictly greater than 
1. One concludes that the derivative at the fixed point ro > 0 is strictly smaller than 1. 
Since this derivative is precisely equal to 2Kd x ^>o(2Kro), this shows that <9 r g(ro,0) < 0. 
So the implicit function Theorem implies that 5 H > r$ is C°°. Using (I1.10|) . one obtains 
that 

r'(0) = d s *sls=o(2JCr 0 ) + r'(0)2Kd x V o (2Kr o ). (D.2) 

Since 2Kd x 'i / o(2A'ro) < 1, the proof of Lemma lD.il will be finished once we have proved 
that (95 'L ( j|,5 = o(2A>o) = 0. One has (recall the definition of Zq in Section lB.ll) 


d 5 y s \ s=0 (2Kr 0 ) = ]T A* 


C cos(0)d s S s k ls =o (0, 2Kr 0 ) d9 
Z 0 (2Kr 0 y 

Jq W cos(6)So(0, 2AVq) „ 


k=—d 


Z 0 (2Kr 0 )* 

Some straightforward calculations show that, for all k = — d,... ,d, 9 € T 


d s Z$\ s=0 (2Kr 0 ) . (D.3) 




( p2it r2ir 

9 j e 2Kr ° cos O) du + 27t J e- 2KroCOs(u hlu 


ue~ 2KroCOs{u) du\ (D.4) 


and 


( /*27t p2tt 

2vr j ' g2iCro cos(0) j ' ( 


-2A> 0 cos(u) dudQ 


+ Zq 

r2n 


(2A> 0 ) u (e 2Kr ° cos W - e - 2 ^ 0 cos( u )^ du j ^ 


= 47TW* 


n2i r r 2 tv 

/ e 2A> 0 cos(0) / e ~2Kr 0 cos (u) du d Q . 

Jo Je 


(D.5) 


Since Ylk=-d^ ku}k = 0, one obtains from (ID. 31) . (ID. 41) and (ID. 51) that <95'L < 5| ( 5 = o(2A>o) = 0. 
This concludes the proof of Lemma ID. 11 □ 

We now turn to the proof of Lemma 17.11 

Proof of Lemma \ r i.l\ Obviously, for 6 £ T, 

qm = qo (e)+5dsqs\ s =o(0) + 0(5 2 ), 

where the error 0{6 2 ) does not depend on 0 £ T. The fact that r -/ (0) = 0 iLemma ID. II) 
implies that dsqs |<5=o (9) only depends on the derivatives of Sg and Zg w.r.t. 6, not w.r.t. 
x. Namely, 

a i, _ d s Si\s =0 (9,2Kr 0 ) d s Zl\ 5=0 (2Kr 0 )S^9,2Kr 0 ) 

ai ® l5 -° (9) “ Z„(2Kr„)2 Z„(2A> 0 )‘ ' 

The expansion found in (17.31) is a simple consequence of (ID.4ft . (ID. 51) and the expression 
of Zq in Section IB. II □ 
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Proof of Lemma | 7. 4 In the case <5 = 0, the projection po defined in (12.151) is given by 
Pq{u) = po(u)(doqo ,..., deqo) = Vo{ u )deqo,nd- Since in this case, the operator Lq = A 
defined in (IB. 41) is essentially self-adjoint in H~^ d (Proposition IB . 1 1) . the projection po 
as a natural representation in terms of the scalar product (■, •)_ 11 i qo d associated to the 
norm defined in (lA.2j) . namely 


(deQ0,nd ! u )_! i/ qo d 

Po U = -o- 

\\®9<10,nd\\-l,l/qo,d 

Using the notations of Section El we deduce that 


(D.6) 


\\deqo,nd 


2 

— 1,1/go ,d 



1 - 


47T 2 


and 


( doqo,nd , u )-x t i/ q ^d 



which precisely gives the first order of (17.51) . The validity of (17.51) comes from the definition 
of the projection Pg in (|2.15l) and the fact that L$ is a relatively bounded perturbation of 
order 5 of the operator Lq = A. □ 
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